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There  will  be  a  workshop  banquet  (also  included  in  your  registration  fee)  held 
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PREFACE 


These  notes  are  intended  to  supplement  the  tutorial  lectures  for  the 
•workshop  by  filling  in  a  few  details  left  out  in  the  lectures.  They  are  not  com¬ 
plete,  but  are  simply  what  was  available  at  the  time.  We  expect  to  have  a 
draft  of  a  complete  set  of  notes  for  the  material  presented  in  the  tutorial  by 
the  end  of  the  year. 

Many  people  contributed  to  these  notes,  although  most  are  probably 
reluctant  to  admit  it  They  began  as  a  very  sketchy  set  of  research  notes 
that  I  developed,  and  have  been  somewhat  rewritten  by  myself  and  the  other 
authors.  In  particular,  Cheng -Chih  Chu  wrote  most  of  Section  2.3.3  on  the 
Solution  of  the  Algebraic  Riccati  Equation  and  helped  throughout  Bruce 
Francis  wrote  most  of  Chapters  2.2  and  2.4  on  Stabilization  and  Approxima¬ 
tion.  Pramod  Khargonekar  also  helped  throughout  but  particularly  with  the 
method  of  proof  used  in  Section  2.3.4  on  Inner-Outer  Factorization.  Gunter 
Stein,  of  course,  was  the  prime  motivator  of  the  whole  tutorial.  Whatever 
seems  enlightening  and  clever  about  the  presentation  is  probably  due  to  the 
other  authors.  Any  obfuscation  is  probably  mine  alone. 

There  are  obviously  lots  of  deficiencies  in  these  notes.  There  are  essen¬ 
tially  no  references,  because  we're  trying  to  take  credit  for  moat  of  the  con¬ 
trol  theory  research  of  the  last  twenty  years.  Actu&Uy,  we  just  didn’t  get  to 
them.  Ditto  on  most  of  the  introductions  to  the  chapters.  The  introductions 
were  supposed  to  motivate  the  technical  details  m  the  subsequent  sections, 
but  hopefully  the  tutorial  lectures  will  do  that.  There  are  probably  lots  of 
typos,  and  whole  sections  that  are  in  the  table  of  contents  haven’t  been  writ- 
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0.0  Notation 


SYMBOL 

USAGE 

• 

*(0:«5*(0 

• 

L  ( x*y)(t )  :*  f  x(t-r)y(r)dT 
— 

* 

2 .  A*  *  complex-conjugate  transpose  of  complex  matrix  A 

3.  P*  *  adjoint  of  operator  T 

u.s 

U(0  s  unit  step  function: 

5(f)  s  unit  impulse 

s 

G(s)  a  two-sided  Laplace  transform  of  g  (t )  j 

1 

| 

ortbogonei  complement  of  X  i 

z 

Z(A)  3  largest  singular  value  of  matrix  A  j 

p 

p(A)  *  spectral  radius  of  matrix  A  I 

Q.L1  Function.  Spaces 


Continuous  Hm»  domain 

Zg(IR.Cm*n):  Hilbert  specs  of  matrix-valued  functions  on  R  ■with  inner 
product 

</.$>:*  7trace[/(t)aj(f}]iS. 

••CO 

$ 

^RC"W):  subspace  of  functions  zero  for  i  <0. 

#a(RCmm>L  :  subspace  of  functions  zero  for  t  >0. 

Pgt  and  Ph£  .  the  orthogonal  projections  from  Z.zORX1*’"*)  onto 
^(RC"m),  ^(RC"***)!  respectively. 

Continuous  frequency  domain  * 

jB:  imaginary  axis. 

Lz(jJR.Cm>m):  Hilbert  space  of  matrix-valued  functions  on  j R  with  inner 
product 

trace^O'y^GO^)]1*"- 

•ffaCyRC"1**):  subspace  of  functions  F{s)  analytic  in  Re  s>0  and  satisfying 
suj  _/* trace  |/r(<y*;y),/,(<7+j6j)](iy  <  •» 
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Lm{j  Wnt  Banach  space  of  (essentially)  bounded,  matrix-v&lusd 

functions,  vdtb  norm 

114=  **  ess^sup  »[p(/«)  ]■ 

^o0S.C“ra):  subspaes  of  functions  F(s)  analytic  &od  bounded  In  Re  * >0. 

/“V, .  tbs  arthsgon&l  projections  tram  lz(jMrtm)sv)  onto 

HaO'E. i?2C?'RCma?,5^  respectively. 


Prefix  B  denotes  real-rational  and  the  prefix  3  dsnotss  the  unit  bail. 
Ilia  symbol  f^***45  denotes  proper  real-rational  matrices.  Sometimes  the 
spaces  are  abbreviated  as  L;(!R),  etc.  or  as  Lz,  etc,  when  context  da  tannines 
the  arguments. 

The  Fourier  transform  yields  the  Mining  [isometric]  isomorphisms; 

4a(R  CW5»)  35  ISC?R  C"“*} 

/Tj(lR.  C***)  a  H^IR,  £w**) 
ffsOR.  cmXi^  hs(j&, 

ii  n 

The  norma  on  these  spaces  are  all  denoted  by  j  j*i  Jy 

A  useful  fact  is  that  the  norm  of  a  matrix  Q  in  Ldj R  Cr‘XT‘)  equa»3  the 
norm  of  the  corresponding  multiplication  operator 

/-*£?/:  C")  **  VJR  C"4); 


that  is. 


G 


S  3UD 


|G/ 


C»).  !!/|Ul 


It  also  equals  the  norm  of  the  operator  restricted  to  H%(j  E,  Cn): 


j | fij !„  k  sup  j  j |<?/  j  ja :  /  £iy8(?E.C%  |j/||^ij 
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0.2.1  Controllability  and  Observability 

Consider  the  system 

x  *  AstBu  .  s(0)  =  0.  (1) 

The  system  or  the  pair  (A.B)  is  eontrollabla  if,  for  each  time  f  t>0  and  final 
state  £i,  there  exists  a  (continuous)  input  it(*)  such  that  the  solution  of  (1) 
satisfies  3  (ft)  =  Sj. 

Theorem  1 

The  following  are  equivalent: 

(i)  (A.B)  is  controllable. 

(ii)  The  matrix  [i?.  AB.  4*5....]  has  independent  rows. 

(iii)  The  matrix  [/4-X/,  B  has  independent  rows  for  all  X  in  C 

(lv)  The  eigenvalues  of  A+BF  can  be  freely  assigned  by  suitable  choice  of  F. 

The  matrix  A  is  said  to  be  stablt  if  all  its  eigenvalues  satisfy  ReX<0.  The 
system,  or  the  pair  (A.B),  is  stabilizabla  if  there  exists  an  F  such  that  A+BF 
is  stable. 

Theorem  2 

The  following  are  equivalent: 

(i)  (A.B)  is  stabilizable. 

(ii)  The  matrix  j<4-XJ\.£  has  independent  rows  for  all  ReXiO. 

We  will  now  consider  the  dual  notions  of  observability  and  detectability 
with  the  system 

x  -  As  ,  x(0)  =  x, 
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y  =  Cx. 


The  -system,  or  the  pair  (C*« 4),  is  obseruoiie  if,  for  every  fi>0,  the  function 
y (f),  t  C  jo.f  xj,  uniquely  determines  the  initial  state  x,. 

Theorem  1‘: 

The  following  are  equivalent: 


(i)  ( C,A )  is  observable. 


(ii)  The  matrix 


C 

CA 

CA2 


has  independent  columns. 


(iii)  The  matrix 


A -XI 
C 


has  independent  columns  for  all  X  in  C; 


(iv)  The  eigenvalues  of  A+HC  can  be  freely  assigned  by  suitable  choice  of  H. 

(v)  (i4',C)  is  controllable. 


The  system,  or  the  pair  (C,A),  is  datactabl*  if  A+HC  is  stable  for  some 
H. 

Theorem  2": 

The  following  are  equivalent: 

(i)  ( C,A )  is  detectable 


/IS\  TV.- 

' 


IA-X/S 

c 


jawf  CCiV 


sna  *11  19* 


(iii)  (A',C)  is  stabilizable. 
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0.2.2  Tranter  Functions 


Consider  the  linear,  time-invariant,  ordinary  differential  equation 
described  by 

x  *Ax  +  Bxl  (1) 

y  *  Cx  ♦  Du 

where  x(t)  €  E  B  is  the  state,  u(f )  e  E  m  is  the  input,  and  y(f )  eE  p  is  the 
output  The  A.B.C,  and  D  are  appropriately  dimensioned  ree\l  matrices. 

Associated  -with  (1)  is  the  convolution  equation 

y(t)*(y*u)(0 

git)  a  Ct*BU(t)  +  DS(t)  (2) 

and,  upon  taking  Laplace  transforms,  the  resulting  transfer  function  is 
y(s)  x  G(s)u(s) 

G{s)  x  C(sI-A)-lB  +  D  (3) 


To  expedite  calculations  involving  transfer  functions  the  notation 


d! 

B 

c 

~D 

i  C(sl-A)~l3  *  D 


(4) 


•will  be  adopted.  Note  that 


[ A  B I 
C  D\ 


is  a  real  block  matrix,  not  a  transfer  func¬ 


tion.  The  product  of  two  transfer  functions  is,  of  course,  the  cascade  of  the 
two  systems  or  just  the  multiplication  of  two  rational  matrices.  The  conven¬ 
tion  will  be  adopted  that  the  product  of  a  matrix  and  a  transfer  function  is  a 
transfer  function  defined  as 


Xu  A_B_  4  XuA  +  XuC\X\\B  +XuD 
Xa  CD  ■  XuA  +  XnCiXnB+XaD- 


A  similar  convention  holds  for  right  multiplication  by  a  matrix 

Suppose  G{s)  is  a  real-rational  transfer  matrix  which  is  proper,  La.,  ana¬ 
lytic  at  s  a«.  Then  there  exists  a  state-space  model  ( A,B,C,D )  such  that 


The  quadruple  ( A,B,C,D )  is  called  a  realisation  of  G.  A  realization  is 
minimal  If  A  has  minimal  dimension.  It  is  a  fact  that  a  realization  is  minimal 
if  and  only  if  ( A,B )  is  controllable  end  (C.A)  is  observable. 


A  basic  object  of  study  will  be  the  transfer  function  and  it  will  be 
assumed  to  have  a  realization.  The  next  section  describes  standard  opera¬ 
tions  on  linear  systems  in  terms  of  transfer  functions  and  their  realizations. 


0.2.3  Operations  on  linear  Systems 


L  Cascade 


Ci 


Bx 

Ci 

By 

At 

Bz 

• 

Cz 

B* 

di! 

Bx 

At 

B2 

Cx\ 

By 

c 2 

Bz 

B\Cz 

B\Dz 

Az  0 

Bz 

* 

0  At 

Bz 

S 

B\Cz  A\ 

B'Bz 

C\  D\Cz 

D\Dz 

t 

D\Cz  C\ 

D\Dz 

Note:  This  realization  may  not  be  minimal. 


2.  Change  of  Variables 

S  -*  S  a  Tx 
y  -»  $  *  Ry 
u  <L  -  Pu. 


A 

B 

A 

A 

T  0 

A 

B 

o 

7 

,  ^ 

C 

\B\ 

TT* 

C 

D 

ii 

0  R 

c 

~D 

1°  Bj 

TAT 1  TBP 

rct1\rdp, 


3.  Stats  Feedback 


u  -*  Q  +  Fx 


A+BF 

C+DF 


B 

PI 

D 

B 


C\D 


I  0 
IF  I\ 


B\ 


4  Output  Injection 

x  a  Ax  +  Bu  -*  £  *  AS  +  Bu  +  Hy 


A 

B 

I  H 

A  B 

A+HC 

B+HD 

C 

~D 

0  / 

C  D 

C 

D 

3.  Tranapoee  (Dual) 

G  -*  Gr 


A 

B 

Ar 

CT 

C 

D 

BT 

dt\ 

6.  Conjugate 


G  -  C? 


A 

B 

-A* 

- cr 

C 

D 

Bt 

Br\ 

7.  Inversion 


Suppose  Dr  is  a  right  (left)  inverse  of  D.  Then  &  = 


A-BDrC 

BDr 

DTC 

& 

a  right  (left)  inverse  of  G. 
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Proof:  The  right  inverse  case  will  be  proven  and  the  left  inverse  case 

follows  by  duality.  Suppose  DDr  -  I.  Then 


G&  * 

A  BETC 

0  A-BETC 

BET 

-BIT 

C  DITC 

DIT 

A  BITC 

BET 

36 

0  A-BITC 

- BIT 

p  ~C~~ 

I 

on  the  right 

yields 


A  0  . 

0 

ccr  a 

0  A-BITC 

•BET 

p  0 

I 

*  / 


Conjugating  the  state 


7  A 

*|o 


on  the  left  and 


i  rl  i -i 

0  /  *  0  I 


Corollary  7  Suppose  Dr  is  a  right  inverse  for  D  and  let 

\a-bdtc\~bz 
0  ~  DrC  I  Z  ■ 

Then 


Gd  ~  DZ. 


Cordlary  T  Suppose  IT  is  »  left  inverse  for  D  and  let 


A-BDrC\-BDr 
~~ZC  j  Z 


Then 


4 


5i?  =  ZD. 


The  following  lemma  characterizes  the  relationship  between  zeros  of  a 
transfer  function  and  poles  of  its  Inverse. 

8.  lemma  Suppose  G  » 

(st,  x,)  such  that 

(X  *  BD  *  5|Z|  , 


X 

5 

c 

13 

with  D  nonsingular.  Then  there  exists 


iff  there  exists  u,  i*0  such  that 


Proof 


G(%) *  0 


(if) 

G(s«}u,  s  0  implies  that  G~l(s)  has  a  pole  at  Thus  3(s,,  x,)  such 
that  Gs,  »*0  and 

(A-BD~lC)x,  *  s,x, 

(ardyif ) 

Set  .  u,  *  -D^Cx,  i*0.  Then 

G(s0  )u,  a  C{stl’-A)~lBut  +  Du^  *  Cx,  -Cr,  *  0 


QED 
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0.2.4  Tinaar  Matrix  Equations  : 

Property  1 :  (Solution  of  Sylvester  Equations)  .  . 

Consider  the  Sylvester  equation 

AX  +  XB  *  C  (1) 

where  A  €  Kru<n,  B  €  3Rw*m,  C  e  K1***1  are  given  matrices. 

Then,  there  exists  a  unique  solution  X  €  if  and  only  If 

*  0.  Vi  *  and  j  *  1 . m. 

Remark : 

In  particular,  if  B  ~  Ar  ,  (1)  is  called  the  "Lyapunov  Equation"  and  the 
nocessary  and  sufficient  condition  for  the  existence  of  unique  solution  will  be 
that  R2[Ai(A)  +  »*  0,  V  ij  a  l, . n. 

Property  2 :  (Solution  of  linear  Equations) 

Consider  the  Unear  equation 

AX*B 

where  A  €  IR71**,  B  e  JR"*"*  are  given  matrices. 

The  following  statements  are  equivalent : 

(i)  there  exists  a  solution  X  e  JR'1*". 

(ii)  the  columns  of  B  €  Rang  a  (A) . 

.  r .  _1  .  f  .  1 

(iu;  ran*  ^  A  if  j  =  ww  |  a  j . 

(iv)  Kar(>ir)CKer(5r). 
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0.3. 1  Gram!  ana  and  Balanced  Realizations 


Suppose 


where 


controllability  grtxmian  Y  as 


A 


is  stable. 


Define  the 


m 

Y  i  feM££,eA'tdt 
and  the  observability  yramian  as 


X  k  feA  t  CCsMdt. 

By  considering  the  corresponding  matrix  differential  equations  it  is 
easily  shown  that  Y  and  X  satisfy  the  Lyapunov  equations 

AY  +  YA'  +  ££'  =  0 
A'X  -  XA  +  CC  *  0 

Note  that  75aO  and  A>0.  Furthermore,  the  pair  (A,£)  is  controllable  ill  Y>0 
and  ( C,A )  is  observable  iff  X>0. 


Suppose  the  state  is  transformed  by  nonsingular  T  to  2-Tz  to  yield  the 
realization 


_  £ 

i  $ 

+ 

A 

E  2\ 

TAT1 

Tff 

G=fc 

)  t) 

with  £ 

?  J&l* 

CTX 

D 

Then  the  g randans  transform  as  ¥  ~  TYT'  and  %  -  (Tl)'XTl.  Note  that 
fx  -  TYXTX  so  the  eigenvalues  of  the  product  of  the  graraians  are  invariant 
under  state  transformation. 

Consider  the  similarity  transformation  F  which  gives  the  eigenvector 


decomposition 
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YX  =  TAT~l  ,  A  =  diag^j....^) 

Then  columns  of  T  are  (possibly  nonunique)  eigenvectors  of  KX"  correspond¬ 
ing  to  the  eigenvalues  {X* }.  It  is  shown  in  Lemma  1  at  the  end  of  this  section 
that  YX  has  real  diagonal  Jordan  form  and  that  A^Q.  This  is  a  consequence  of 
YstQ  and  AfeO. 


Although  the  eigenvectors  are  not  unique,  In  the  case  of  a  minimal  reali¬ 
zation  they  can  always  be  chosen  such  that 

?-  777'  =  £, 

X  -  cr'yxr1  *  E. 


where  £  *  diag(ffltffa . crn)  and  £8  «  A.  This  new  realization  will  be 

referred  to  as  a  balanced  realization  (also  called  internally  balanced). 


Suppose  G  s 


A 

B 

C 

~D 

is  a  balanced  realization  for  G  and  can  be  parti¬ 


tioned  as 


Au  A 12 

Bx 

Aai  Ass 

Bz 

Cx  C2 

D 

with  corresponding  partitioning  of  the  balanced  grumian  E 


o; 

0  Eci 


.  Sup¬ 


pose  El=diag(ffj,ffa,...ffr).  £s-diag(jtM.l,dr+a . °w)  sad 


ffxJ&ffzS:  •  •  5sar>ffr^la:iTt.+8^  •  •  •  Then  it  is  immediate  that  the  truncated 


system 


r. 

Km 

O  s 

pn 

D 

i3  balanced  since 
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^11^*1  ^i^n'  s  0 

Axi'Di  *  ^l^u  ■**  ^i'Ci  *  0* 

It  can  also  be  shown  (Silverman  and  Pemebo)  that  &  minimal  realisation  for 
Gr  is  stable,  although  in  certain  (non-generic)  cases  Au  may  have  uncontroll¬ 
able  or  unobservable  ;v-axis  eigenvalues. 

Lemma  1  Product  of  Positive  Semi-Definite  Matrices  is  Similar  to  a  Positive 
Semi-Definite  Matrix 

Proof:  Let  X  and  Y  be  positive  semi-definite.  First  perform  an  orthogonal 
transformation  so  that 

A,  o  1  Yu  Ytt ' 

*  -*  [o  oj  A'>0  dia*oaa1’  Y  [ru- 

AiFji  Aj Yia 

By  this  transformation  XY  is  similar  to  q  q  Now 

Ax Yn  At Ya  A?  0  ]( A^uA?  A^u  |fAf*  0 

,0  0  *  o  /][  o  o  Jl  0  /]• 

and  it  is  easy  to  find  a  matrix  Z  such  that 

ApruAl*  0 1  | /  z\k$YnAt  $YU \\l  -Z 

0  oj  *  [o  /)[  0  0  j[o  / 

(Z  exists  because  the  columns  of  A^Fjj  span  the  columns  of  A^Fu  owing  to 
the  fact  that  Y  is  positive  semi-definite).  The  left  band  side  of  this  last  equa¬ 
tion  is  positive  semi-definite  and  similar  to  XY.  If  ?CY*Q  it  is  possible  to  find 
a  matrix  Tt  such  that 
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Tr'ABT ,  =  A  =  diasp, . Xjr.0 . 0)J 

where  AjSsAgat  •  •  •  SeAjr>0. 

Q.E.D. 

Now  consider  two  gramians  X  and  Y.  Let  us  suppose  XY*§,  so  that  T\ 
can  be  chosen  as  above: 

TrlXYTi  a  A  =  diagf(X, . Ajr.O . 0)1 


•where  Aja^t  •  •  •  SeXj^’O.  Under  this  transformation  the  gramians  become 

q  *  rrWrl)\  r  *  jyjtx. 

and  QR-L  Because  Q  and  R  are  symmetric,  RQ=A‘=A=QR  and  so  Q,R  and  A 
commute.  Both  Q  and  R  must  leave  the  eigenspaces  of  A  invariant  and  so 
are  of  the  form 


Q  »  diag  \Q , . 3 


E 


R 


where  3  is  a  square  matrix  whose  size  equals  the  dimension  of  the  A^  eigen* 
space  of  A  and  EF= 0  where  E  and  F  are  square  matrices  the  size  of  the  ker¬ 
nel  of  A  Of  course,  ail  the  3's  are  summetric  so  it  is  possible  to  find  an 
orthogonal  matrix 

V  =  diao  Iw . tt.M*,! 

'  'J 


such  that  W~lQ(1Y~1)'  and  are  diagonal.  Note  that  this  same  W 

gives  a  diagonal  W'RW  and  leaves  A  alone. 


5 

The  transformation  Ti-T^W  diagonalizes  both  grami&ns.  It  is  now  obvi* 
ous  how  to  construct  Tj  so  that  Ts  'X(T «T1)’  and  7' j' YT j  are  diagonal  and  the 
controllable  and  observable  portions  are  equal. 
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0.3.2  Inner  Transfer  Functions 


Let  G  — 


dl 

B 

c 

D 

.  Then  G  La  inner  if  GmG~I  end  co-innar  if  GG*-I.  Note 


that  G  need  not  be  aquare.  Inner  end  co-inner  ere  dual  notions  and  are  often 
called  all-paaa. 


If  G  6  ft£xm ,  p>m  la  inner  then  any  ia  called  a  comple¬ 

mentary  inner  f&Otor  (CIF)  if  Crjj  is  square  and  inner.  The  dual  notion  of 
complementary  co- inner  factor  is  defined  in  the  obvious  way. 

The  following  lemma  is  useful  in  characterizing  inner  transfer  functions 
in  terms  of  a  realization 

Lemma  1.  Suppose  2  X=X  e  R**"  such  that 

i)  A'X  *  XA  *  CC  -  0 

ii)  B’X  +  D‘C  =  Q 

Then  GaG  —  D'D. 

Proof:  Suppose  that  i)  and  ii)  hold.  Then  conjugating  the  state  of 


A  0 

B 

GmG  = 

-C'C  •A' 

- CD 

D'C  B' 

D'D 

hy 


I  0 
-X  I 


on  the  left  and 


I  0 

-l 

/  0 

-X  I 

X  I 

on  the  right  yields 


A  0 

B 

GmG  = 

-(A'X+XA+CC)  •A' 

< XB+CD ) 

B'X+D'C  B' 

D'D 

Now,  applying  i)  and  ii)  yields 
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A  0 

B 

G*G  = 

0  -A' 

0 

[O  B' 

D’M 

a  D'D 


By  duality,  we  have  the  following 

Lemma  1*  Suppose  3  T  e  JR"*"  such  that 

i)  AY  +  YAt  +  BB'  -  0 

ii)  CY  +  DB'  *  0 

Then  CG*=DD‘. 


These  two  lemmas  immediately  lead  to  one  characterization  of  inner 
matrices  in  terms  of  their  state  space  representation.  Simply  add  the  condi¬ 
tion  that  D'D-I  (DD'=I)  to  lemma  1  (1’)  to  get  GmG-I  {GG*-I).  Further¬ 
more.  by  adding  a  few  additional  assumptions,  the  conditions  in  the  lemmas 
become  necessary  as  well  as  sufficient.  This  leads  to  the  following  complete 
characterization  of  stable  inner  transfer  functions  in  terms  of  a  minimal 
realization. 


Suppose  G- 


A  B\ 
C  D\ 


is  stable  and  minimal. 


Then  the  gramians  X  and  Y 


satisfying 


A'X  +  XA  w  CC  -  u  (2) 

AY  +  YA‘  +  3B’*0  (3) 


exi«t  and  are  unique. 


Corollary  1  G  is  inner  ill 
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i)  B'X  +  PCs  0 

ii)  D’D  =  / 

Corollary  1'  G  is  co-inner  iff 

i)  CY  +  DB' s  o 
U)  DP  =  / 


Proof  Sufficiency  of  i)  and  ii)  follows  immediately  from  the  lemmas.  For 
necessity,  suppose  £?•£?*/.  From  1)  and  2)  this  implies  that 


'  A 

B 

B'X+PC 

0 

(4) 


PD  =  / 


(5) 


Since  (A, 5)  is  controllable,  (4)  implies  that  B’X+D'C=Q.  The  co-inner 
case  follows  by  duality. 

This  characterization  of  inner  transfer  functions  plays  a  central  role  in 
the  synthesis  theory.  It  allows  the  construction  of  inner  transfer  functions 
by  solving  algebraic  equations. 
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0.4.1  linear  Fractional  Transformations 


Suppose  P  =  |£“  ,  ic^.jre  We 

will  adopt  the  notation 

Ti(P,K)  4  Pn  +  Pi2K(I-PaKj-lPn  (1) 

and 

Tu(PA)  i  Pzt  +  PiAV-PiAr'Pu  (2) 


The  linear  fractional  transformations  (LFT)  are  illustrated  in  Figure  1. 
The  i  denotes  that  the  second  argument  is  fed  back  in  the  lowtr  block,  and 
the  u  denotes  feedback  in  the  upper  block. 


An  important  property  of  LFTs  is  that  any  interconnection  of  LPTs  is 


again  an  LFT.  Suppose  J 


T i(P'Ti(J.Q))  =  W,Q)  (3) 

Tu(j,ru(p,D)  =  fu(t.a)  (4) 


where 


T  = 


Tn  Tiz 
Tzi  Tn. 


Pm  *  P viJ ii(I~PzzJ ii)">1^52:  P \zU“J n Ptt)~lJ iz 

J2l{I~PzzJ  ll)~XPz\  Jzi  +  JziPz2{I~JllPzd~'JlZ, 


(5) 


Equations  (3)  and  (4)  are  illustrated  in  Figure  2.  Note  that  if  Ti{J,Q)  is  a 
parametrization  of  a  controller,  Ft (T.Q)  is  affine  if  and  only  if  Tz 2=0.  This 
type  of  controller  parametrization  will  play  an  important  role  in  the  syn¬ 
thesis  theory. 


I^+/o  du  cf-1* 


This  part  of  the  notes  describes  recent  results  on  the  problem  of  analyz¬ 
ing  the  performance  and  robustness  properties  of  systems.  We  believe  that 
tue  approach  described  is  providing  the  foundation  for  a  new  paradigm  for 
catrol  theory  broader  in  scope  and  content  than  that  of  Classical  or  Modem 
.tmtrol  Theory.  An  important  aspect  of  this  new  paradigm  is  the  treatment 
it  -jives  to  model  uncertainty. 

•  Modem  Control  Theory,  the  dominant  paradigm  for  the  past  20  years, 
u&s  its  basis  in  Stochastic  Optimal  Control  and  Estimation  Theory.  This 
theory  essentially  restricts  model  uncertainty  to  additive  noise.  The  theory 
provides  a  methodology  for  analyzing  the  impact  of  noise  on  system  perfor¬ 
mance  and  synthesizing  to  reduce  that  impact 

The  inadequacies  of  this  view  of  uncertainty  became  widely  accepted  in 
the  Lite  1970’s,  as  robustness  to  plant  uncertainty  became  a  major  theme  in 
the  Modem  Control  Theory  community.  Ironically,  this  involved  a  renewed 
interest  in  the  Classical  Control  paradigm  which  Modem  Control  displaced 
within  the  theoretical  community  (if  not  among  practicing  engineers).  This 
new  direction  provided  useful  design  tools,  Including  Singular  Value  Analysis 
and  Multivariable  Loop  Shaping. 

While  providing  an  important  perspective,  as  well  as  practical  tech¬ 
niques,  the  methods  based  on  singular  values  still  require  rather  restrictive 
assumptions  about  uncertainty.  In  particular,  plant  uncertainty  must  essen¬ 
tially  be  modelled  as  a  single  "unstructured  perturbation." 

The  Structured  Singular  Value  (5SV),  (jl,  was  developed  several  years  ago 
to  correct  this  deficiency  in  singular  values.  In  the  context  of  the  general 
framework  discussed  in  this  memo,  the  SSV  provides  a  very  powerful 
mathematical  tool  for  the  analysis  of  complex  systems.  Indeed,  we  believe 
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that  this  framework  together  with  the  SSV  and  the  synthesis  techniques  dis¬ 
cussed  later,  has  the  potential  to  form  the  basis  for  a  new  paradigm  for  con¬ 
trol  theory. 

The  remainder  of  this  part  of  the  notes  describes  the  general  framework 
for  control  system  analysis  and  synthesis  which  includes  all  the  viewpoints 
discussed  as  special  cases,  in  particular,  the  assumptions  about  uncertainty 
required  by  each  methodology  are  compared.  In  this  context,  the  words 
analysis  and  synthesis  have  specific  meanings. 

Analysis  is  used  to  describe  the  process  of  determining  whether  a  given 
system  has  the  desired  characteristics.  In  general  this  may  range  from  the 
use  of  mathematical  tools  to  simulation  to  experimentation,  although 
analysis  is  typically  applied  primarily  to  describe  the  former.  Synthesis,  on 
the  other  hand,  is  the  process  of  finding  a  particular  system  component  to 
achieve  desired  characteristics,  which  are  typically  expressed  in  terms  of 
some  analysis  tools.  Analysis  and  synthesis  are  just  two  aspects  of  the  more 
general  problem  of  engineering  design. 

The  discussion  which  follows  first  considers  analysis,  then  briefly 
touches  on  synthesis  and  ends  with  some  illustrative  examples.  The  next 
part  on  Synthesis  Theory  will  take  up  that  question  in  more  detail. 
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1.2.1  General  Framework 

Various  modelling  assumptions,  will  be  considered  and  the  impact  of 
these  assumptions  on  analysis  and  synthesis  methods  will  be  explored.  Con¬ 
sider  the  diagram  in  Figure  L  This  is  the  general  framework  to  be  con¬ 
sidered.  Models  of  this  form  are  typically  constructed  from  components 
which  also  have  this  form.  The  nominal  model  provides  the  basic  intercon¬ 
nection  structure  between  the  signals,  perturbations  and  controller,  as 
shown.  It  has  three  inputs  and  outputs,  each  consisting  of  a  vector  of  signals. 

As  typical  examples,  consider  the  following  filtering  and  control  prob¬ 
lems.  First,  a  simple  filtering  problem  is  given  in  the  diagram  in  Figure  2. 
This  may  be  rearranged  as  shown  in  Figure  3  to  fit  the  general  framework.  In 
order  to  simplify  the  diagram,  no  perturbation  was  included. 

A  typical  control  problem  might  look  like  the  diagram  in  Figure  4  where 
again,  for  simplicity,  no  perturbations  are  included.  This  too  can  be  rear¬ 
ranged  to  fit  the  general  framework,  although  the  diagram  is  complicated. 

Any  system  may  be  rearranged  to  fit  the  form  of  this  general  frame¬ 
work.  Although  the  interconnection  structure  can  become  quite  complicated 
for  complex  systems,  many  software  packages  are  available  which  could  be 
used  to  generate  the  interconnection  structure  from  system  components. 

Note  that  uncertainty  may  be  modelled  in  two  ways,  either  as  external 
inputs  or  as  perturbations  to  the  nominal  nodeL  The  performance  of  a  sys¬ 
tem  is  measured  in  terms  of  the  behavior  of  the  outputs  or  errors.  The 
assumptions  which  characterize  the  uncertainty,  performance  end  nominal 
model  determine  the  analysis  techniques  which  must  be  used. 

The  most  fundamental  assumption  that  is  made  throughout  is  that  the 
nominal  model  i3  a  finite  dimensional  ordinary  differential  equation  and  is 
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linear  and  time  invariant  (LTIODE).  The  uncertain  inputs  are  assumed  to  be 
either  filtered  white  noise  or  weighted  Lp  signals.  Performance  is  measured 
as  weighted  output  variance  or  weighted  output  Lp  norm.  The  perturbations 
are  assumed  to  be  themselves  LTIODE's  which  are  norm-bounded  as  input- 
output  operators.  Various  combinations  of  these  assumptions  form  the  basis 
for  all'the  standard  linear  systems  analysis  tools. 

Given  that  the  nominal  model  is  an  LTIODE,  the  interconnection  system 
has  the  form 


Pn 

P 12  Pis 

P  a 

P21 

,  P32  Pa 

P* 

Pfu  Pa 

A 

B  i  Ba 

b3 

Ci 

Du  D a 

Da 

*• 

Ca 

D%\ 

Da 

c. 

Dsi  Da 

Da 

and  the  total  system  is  a  linear  fractional  transformation  on  the  perturbation 
and  the  controller  given  by  * 


e  =  Tu(Fi(P,K),L)  u 

=  (2) 
Since  the  focus  of  the  current  discussion  is  on  analysis  methods,  the 
controller  may  be  viewed  as  just  another  system  component  and  absorbed 
into  the  interconnection  structure.  Thus  the  analysis  framework  reduces  to 
the  diagram  in  Figure  5  where 


a  =  Tu(PA)u 

~  |^22  +  Pz\&(I ~P \1&)~^P  12]  ^ 


(3) 
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Note  that  the  P’s  in  (2)  end  (3)  are  not  necessarily  the  same.  Table  1 
and  the  discussion  which  follows  summarize  the  various  assumptions  and 
resulting  analysis  and  synthesis  tools.  In  each  case,  stability  of  the  nominal 
must  be  evaluated.  Since  P  is  assumed  to  have  the  state-space  representa¬ 
tion 


A 

8 1  Bz 

p  « 

C , 

Dn 

Ca 

Du  Dj2 

this  may  be  done  by  checking  that  all  eigenvalues  of  A  lie  in  the  open  lhp. 
There  are  alternatives  to  this  approach  but,  for  simplicity,  it  will  be  assumed 
that  the  nominal  plant,  with  controller  is  closed  loop  stable  in  the  sense  that 
all  eigenvalues  of  A  are  in  the  open  lhp. 

Given  nominal  stability,  the  entries  in  the  table  may  be  interpreted  as 
filling  in  the  following  general  performance/robustness  theorem; 


General  Analysis  Theorem  (GAT) 

Given 


Input  Assumptions 


and 


Than 


Perturbation  Assumptions 


Performance  Specification 


if  and  only  if 


Analysis  Test 


The  details  of  each  case  ■will  be  considered  in  the  following  sections. 


> 


rn&a-Surtd 

Oicfpa+S 


£.omnrtan<le.dt 
V<xr»  curies 


fncasurem&At 

noise. 
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1.2.2.  Stochastic 


Case  la  involves  unit  covariance  white  noise  input  with  output  variance 
as  the  evaluation  criteria.  Since  no  perturbation  is  allowed,  the  problem 
reduces  to  the  diagram  in  Figure  1  and  Ft#2  s)  *  Fjs  .  Note  that  colored 

noise  or  weighted  variance  could  be  used  as  shown  in  Figure  2.  This  reduces 

to  the  general  case  by  absorbing  the  weights  and  into  Pss.  as 

Pn  -  WzGVv  In  practice,  it  is  essential  to  use  weights  to  reflect  spatial  and 

* 

frequency  variations  in  inputs,  perturbations  and  output  specifications,  but 
in  every  case,  these  weights  may  be  absorbed  into  nominal  modeL 


In  Case  lb  the  input  is  an  uncertain  delta  function,  which  is  equivalent 
to  uncertain  initial  conditions.  The  performance  specification  is  the 
expected  value  of  the  Lz-norm  of  the  output 

Case  1  forms  the  foundation  of  Stochastic  Optimal  Control  Theory.  Case 
la  includes  the  standard  linear  stochastic  filtering  problem  and  Case  lb 
includes  the  standard  linear  quadratic  optimal  control  problem.  These  ere 
combined  to  obtain  the  full  LQG  problem,  which  is  again  Case  la.  These 
assumptions  and  resulting  analysis  methods  have  been  the  dominant  para* 
digm  in  the  control  community  for  over  20  /ears. 

The  development  of  this  paradigm  has  stimulated  extensive  research 
efforts  and  been  retponsible  for  important  technological  innovation,  particu¬ 
larly  in  the  area  of  estimation.  The  theoretical  contributions  include  a 
deeper  understanding  of  linear  systems  and  improved  computational 
methods  for  complex  systems  through  state-space  techniques.  The  major 
limitation  of  this  tu»ory  is  the  lack  of  formal  treatment  of  uncertainty  in  the 
plant  itself.  By  allowing  only  additive  noise  for  uncertainty,  the  stochastic 
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theory  ignored  this  important  practical  issue.  Plant  uncertainty  is  particu* 
l&riy  critical  in  feedback  systems. 
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1.2.3  hm  Frequency  Ennuis.  Methods 


Case  2  involves  an  attempt  to  correct  some  of  the  deficiencies  of  Case  1 
by  moving  to  an  unknown  but  bounded  (in  an  £g  sense)  framework.  This 
allows  both  types  of  plant  uncertainty  to  be  handled  in  a  common  framework, 
albeit  in  a  limited  manner. 


Case  2a  is  an  ig  version  of  Case  la.  The  input  is  constrained  to  lie  in 
Biz  as  a  time  signal  (unit  ball  in  Lg)  and  the  performance  is  specified  in 
terms  of  the  output's  L g  norm.  With  no  perturbation,  the  analysis  test 
involves  simply  the  L%  induced  operator  norm.  Le.  Lm  on  the  transfer  func¬ 
tion  Pa. 

The  CAT  in  this  case  is 


Theorem  1 


iff 


2  S  1  for  all 


*  1 


Pj2 


!  £  1 


Although  this  theorem  is  a  trivial  restatement  of  the  definition  of  induced 
norm,  it  means  that  th e-analysis  test  is  on  exact  characterization  of  the  per¬ 
formance  requirement. 

Case  2b  is  significant  departure  from  the  previous  three.  It  involves 
maintenance  of  stability  in  the  presence  of  perturbations.  The  block 
diagram  for  7U'(P,A)  is  shown  in  Figure  1.  There  are  many  ways  to  state  the 
GAT  for  this  case,  depending  on  the  desired  notion  of  stability  and  assump¬ 
tions  on  A.  The  distinctions  are  somewhat  subtle,  but  are  important  from  a 
theoretical  point  of  view.  Nevertheless,  they  do  not  significantly  imp  act  the 
application  of  the  theory. 
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The  A‘s  are  assumed  to  be  LTIQDE's,  so  that  A  e  RHm-  The  assumptions 
A  €  Cor  A  £  C Hm,  give  the  same  result  The  distributed  case,  A  e  Hm,  causes 
some  additional  technical  difficulties  and  is  not  the  focus  of  these  notes. 

The  fo&owing  theorem  treats  internal  stability 


Theorem  2  TU'(P,A) 


is  internally  stable  for  all  A  e  BRHm 


Note  that  input-output  stability  of  FU(P,A)  is  not  necessarily  the  same  as 
internal  stability.  In  particular,  the  following  statement  is  not  true: 


Not-A- Theorem 


2<»  for  all 


stfil  and  A  e  BRH. 


iff 


11 


ts  1 


Counterexample  Suppose  I  ML  >1  but  Pz\i 


sO. 


From  now  on  stability  will  mean  infernal  stability,  but  be  denoted  by 
2  <  «  in  the  table,  even  though  this  is  definitely  an  abuse  of  notation. 
Note  that  genetically  this  distinction  between  internal  and  1-0  stability  does 
not  exist. 


As  in  Case  1.  it  is  essential  to  allow  weights  on  inputs,  outputs  and  per¬ 
turbations.  As  before,  these  weights  may  be  absorbed  into  the  nominal 
model.  This  allows,  without  the  loss  of  generality,  the  use  of  signals  and  per¬ 
turbations  which  are  in  unweighted  unit  balls.  Thus  implementation  of  the 
analysis  tools  requires  oni,  a  method  for  constructing  interconnected  sys¬ 
tems  and  a  method  for  evaluating  the  appropriate  norm.  The  former  applies 
to  all  cases,  whereas  the  latter  requires  a  different  norm  in  each  case. 
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Note  that  in  Case  2  both  uncertain  inputs  and  uncertain  plants  can  be 
handled  with  the  same  analysis  tool.  This  approach  is  particularly  useful  for 
feedback  problems  where  both  types  of  uncertainty  have  significant  impacts 
on  system  performance.  Case  2  has  attracted  a  great  deal  of  research 
interest  recently,  and.  is  currently  the  popular  new  paradigm  in  the  mul¬ 
tivariable  control  community.  Although  implicit  in  the  methods  of  classical 
control  and  some  more  modem  work  (e.g.  Zanies'  conic  sector  theory  circa 
early  60's  and  Horowitz's  80’s  work),  the  approach  did  not  gain  wide  atten¬ 
tion  until  the  late  '70's. 

The  current  interpretation  is  a  consequence  of  research  done  in  the  late 
'70' s.  (Doyle  and  Stein,  Safonov  etc).  This  interpretation  involves  singular 
values  as  an  analysis  method  and  singular  value  loop  shaping  as  a  synthesis 
approach.  The  so-called  LQG  Loop  Transfer  Recovery  (LQG/LTL;  (Stein  and 
Doyle)  combines  the  synthesis  methods  of  Case  1  with  the  analysis  methods 
of  Case  2  to  produce  a  hybrid  synthesis  method.  This  gives  an  ad  hoc 
approach  to  Cast;  2  that  can  be  effective  for  many  multivariable  problems. 

■“  Another  approach  to  synthesis  for  Case  2  is  the  so  called  Hm  or  I-/  Hm 
methods  introduced  to  the  control  community  by  Zames  and  Helton  (with 
additional  contributions  by  Francis,  Pearson,  Glover  etc.).  The  Lm/ Hm 
methods  for  Case  2  are  analogous  to  the  Lz/  Hz  methods  of  Case  1  with  the 
exception  that  for  Case  2  the  Lm  rather  than  Lz  norm  is  optimized.  The  solu¬ 
tion  to  the  general  Lm/  Hm  problem  will  be  presented  in  the  Synthesis  part  of 
these  notes. 

The  main  objection  to  Case  2  is  the  resuictive  assumptions  about  uncer¬ 
tainty  (recall  this  was  also  the  objection  to  Case  1).  Although  case  2  allows 

♦ 

both  uncertain  inputs  and  perturbations,  analysis  can  be  performed  for 
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either  individually  but  not  both  together.  Thus  a  system  can  be  shown  to 
remain  stable  when  perturbed  and  have  acceptable  response  to  uncertain 
inputs  when  A  =  0  but  response  when  A*0  is  not  known.  Only  crude  bounds 
can  be  obtained  with  the  methods  of  Case  2. 

An  additional  limitation  of  Case  2b  is  that  all  plant  uncertainty  must  be 
modelled  as  a  single  norm-bounded  perturbation.  Typically,  uncertainty  is 
present  throughout  a  system.  Suppose  that  a  system  is  built  from  com¬ 
ponents  which  are  themselves  uncertain  and  that  component  uncertainty  is 
modelled  as  norm-bounded  perturbations.  This  situaton  can  tea  rearranged 
to  fit  the  general  framework  but  the  perturbation  for  the  total  system  has 
structure.  The  problem  of  structured  uncertainty  is  taken  up  in  the  next 
chapter. 
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1.3.1  Introduction 


This  chapter  considers  the  problem,  of  stability  with  structured  uncer¬ 
tainty  and  of  performance  in  the  presence  of  structured  uncertainty.  Typi¬ 
cally,  uncertainty  is  present  throughout  a  system.  Suppose  that  a  system  is 
built  from  components  which  are  themselves  uncertain  and  that  component 
uncertainty  is  modelled  as  norm-bounded  perturbations.  This  situaton  can 
be  rearranged  to  fit  the  general  framework  but  the  perturbation  for  the  total 
system  has  structure.  This  can  be  seen  schematically  in  Figure  L 


Note  that  the  interconnection  model  P  can  always  be  chosen  so  that  A  is 
block  diagonal,  and  by  absorbing  any  weights,  |  A  .<1.  The  results  of  Case 
2b  can  be  applied  in  two  ways: 

1)  iiFuji  <1  implies  stability,  but  not  conversely.  This  can  be  arbitrarily 


l!i» 

conservative,  in  that  stable  systems  can  have  arbitrarily  large 


u 


2)  Test  for  each  A(  individually.  This  can  be  arbitrarily  optimistic  because 
it  ignores  interaction  between  the  A(. 


The  difference  between  the  bounds  obtained  in  1)  and  2)  can  be  arbi¬ 
trarily  far  apart.  Only  when  they  are  close  can  conclusions  be  made  about 
the  general  case  with  structured  uncertainty. 

These  two  limitations  of  Case  2  (and  l)  have  motivated  much  of  the 
research  described  in  these  notes.  The  result  is  a  new  paradigm  described  in 
Case  3.  The  problem  in  Case  3  involves  exactly  that  of  structured  uncer¬ 
tainty. 


Consider  the  system  in  Figure  2.  Stability  and  performance  analysis  of 
this  system  requires  a  new  matrix  function,  the  structured  singular  value 
(SSV),  denoted  by  (jl  Before  proceeding  with  Case  3.  a  digression  to  discuss  fx 


2 


will  be  taken.  For  details,  see  the  reprints  which  accompany  this  writeup. 
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1.3.2  SSV for  Constant  Matrices 

The  problem  is  to  test  for  for  sets  of  A.  Two  standard 

results  are 


1)  ctat(/-ifA)*0  Vie  A  p(A)<l 


Iff  U(M)  se  1 


2)  dat(I-ML)*Q  v  A  e 


A/  A  €  C.  X 


<1 


Iff  p(if)  s  1  where  p(M)  3  raaac  jA((Af)  | 

As  a  generalization,  consider  a  function  p.  with  the  properties  that 
p(aJl)=  |a|/i(tf)  and 

3)  dat(I-&Ib.)*Q  v  AS  Lo^(Ai.Ab . A«)  ^(At)*! 


iff  n(M)  <  1 


Obviously,  p.  is  a  function  of  M  which  depends  on  the  structure  of 


A 


To  be  precise,  a  multi-index  could  be  constructed  which  would  specify  use 

[  1 

structure  of  |  A  j  and  p  would  depend  on  that  index  For  this  informal  dis¬ 
cussion,  just  keep  in  mind  thi3  fact  and  assume  that  a  structure  is  specified. 
Clearly  ff  and  p  are  special  cases  of  p  for  particular  structures  as  indicated 
above.  Furtnermore,  tor  any  structure 


p(M)  *  p{M)  S  73{M).  (4) 

Given  these  bounds,  how  important  is  p'i  The  answer  can  be  clearly  seen 


from  the  following  examples: 
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<$!  0 

Suppose  A  ~  q  and  consider 

1)  Si-  q  £  p(tf)=0  V(JI)=1 

dat(I-£iL)-l  so  p,(tf)*0 
_u  u 

2)  p(M)~ 0  ?(JV)=1 

d*i(/-^A)  =  1  +  ~ so 

Thus  neither  p  nor  3  provide  useful  bounds' even  in  simple  cases.  The 
only  time  they  do  provide  reliable  bounds  is  when  p*V.  Thus  better  bounds 
on  p  are  needed  to  pursue  the  problem  in  Case  3. 

•  For  the  rest  of  the  discussion  fix  a  structure  for  the  A’s  as 


where  the  sets  SI  and  Q  match  the  structure  of  X-  Note  that  the  SI  and  Q. 
leave  X  invariant  in  the  sense  that 
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1)  A  e  X  u  e  U  imply  9(AU)  =  9(  UL)  =  9(A) 

2)  A  eX  DC  R  imply  DLD~X  =  A 

From  these  two  properties  and  the  definition  above  expression  for  fx,  one 
immediately  obtains 

p(tf  0)  a:  i4M)  *  jttf  V(DMJrl) 

The  first  important  theorem  about  /x  is 
Theorem  1  npax  p(MU)  -  pI^M) 

This  theorem  expresses  p.  in  terms  of  familiar  linear  algebraic  objects. 
Unfortunately,  the  implied  optimization  problem  is  nonconvex  so  it  does  not 
immediately  yield  a  computational  approach.  The  second  important 
theorem  is 

Theorem  2  If  rue 3  piM)~  inf^  T3{DMD~X) 

This  theorem  states  that  if  there  are  3  or  fewer  blocks  (  no  restriction 
on  size),  then  tx{M)  is  just  9  of  a  block  diagonal  similarity  of  M.  Furthermore 
V(DMD~V  is  convex  in  D  so  that  the  inflmum  can  be  found  by  search  over 
n-1  real  parameters. 

The  theorem  is  not  true  for  n:s4,  but  it  is  conjectured  that 
^inf^  TS{DMD~X)  still  provides  a  reasonably  tight  bound  for  /i.  Also,  many 

problems  of  interest  have  3  or  fewer  blocks  so  this  provides  a  reasonable 
computational  scheme. 

Another  important  aspect  0*  thii  theorem  IS  that  p,  may  be  vas~5d  as  9 
plus  scaling.  Thus  the  general  synthesis  methods  recently  developed  to 
optimize  the  X.  norm  (i.e.  9)  may  be  applied,  via  scalings,  to  optimize  fx. 
This  will  be  discussed  more  in  the  syhthesis  section.  Now  back  to  Case  3. 


1 


1.3.3  SSV Analysis  of  Systems 


Abuse  notation  and  define 


IHU  *  sup  uilHju)).  (1) 

Although  j*  is  not  a  norm,  this  will  be  convenient  Recall  tate  ||#||^ 
ii  a  function  of  U  which  also  depends  on  the  assumed  structure  of  the  per¬ 
turbations. 


Case  3a  involves  stability  in  the  presence  of  structured  perturbations 
tnd  the  result  is  analogous  with  Case  2b.  In  fact  3a  reduces  to  2b  in  the  case 
that  there  is  a  single  block  in  the  perturbation.  Suppose  that  A  E  BRHm  and 
the  A’s  have  the  structure  A=diap(Ai,Aa,Ao,  ....  An).  The  CAT  for  Case  3a  is 


Theorem  1  YU(P,L)  is  internally  stable  for  all  structured  A  £  BRHm 

m  i>,4  *  i 

« 

Case  3b  puts  everything  together  and  is  really  the  payoff  for  analysis.  The 
)  rohlem  is  to  check  that  jleJusSl  is  satisfied  for  all  u  {£1  and  all  struc¬ 
tured  perturbations.  Recall  that  from  2a  and  2b  that  both  stability  with  a 
single  perturbation  and  performance  with  Lz  inputs  involve  the  same  test 

; :  M 

using  VL  although  on  different  parts  of  the  system.  This  means  that  the 
system  in  Figure  1  has  internal  stability  and  e  j  SI  for  all  u  Z.<1  and 

A=sKay(Aj,Uj . An)  £  BRHm  if  and  only  if  the  system  in  Figure  2  has 

interned  stability  for  all  structured  A  and  all  A» *i  £  BRHm.  Thi3  i3  exactly 

Case  3a  with  the  structure  A=tfiaj(Ai,Aj . An.An+i).  Using  this  structure 

for  u  yields  the  following: 
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Theorem  **-' 

?U(P,&)  i a  tofeamally  stable  and 
A— cfici^  (Aj.Ag,  .  . .  .A*)  £ BRHm 


:  aSl 


for 


i£  PjttsSl 


and 


This  is  a  ramark&bly  useful  theorem.  It  says  that  ||p|j^ss  1  implies  not 
joly  stability  for  all  structured  perturbations  but  also  that  |  jfl  1 12  5  1  for  all 


j'u|[as£l  and  ail  structured  perturbations.  Furthermore.  |  J/3) ^  >  1  implies 
that  there  exists  a  u  with  |  |uj  js  :£  1  and  a  structured  A  such  that  either 
jj«j|a  >  1  or  FU(P,A)  is  internally  unstable.  This  is  the  first  general  result 
which  guarantees  performance  for  a  whole  sat  of  plants  and  gives  an  exact 
(nonconservative)  analysis  test 
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1.4.1  A  Glimpse  at  Synthesis 

This  will  be  a  sketchy  outline  of  the  new  synthesis  results.  The  details 
are  somewhat  complicated  and  are  treated  in  Pert  2  which  is  devoted  to  the 
synthesis  theory.  At  this  point,  we  simply  want  to  point  out  how  the  analysis 
theory  discussed  in  this  part  leads  naturally  to  certain  synthesis  questions. 

From  the  analysis  results,  we  see  that  each  case  boils  down  to  evaluating 

\\Pij\l  «=2.«  or  n  (1) 

for  some  transfer  function  Pq.  Thus  when  the  controller  is  put  back  into  the 
problem,  it  involves  just  a  simple  linear  fraction  transformation  as  shown  in 
the  diagram  in  Figure  1.  (Note:  the  Py's  here  are  not  the  same  as  the  Py’s 
in  the  previous  sections) 

Each  case  then  leads  to  the  synthesis  problem 

nun'iFiCPJo!1  for  a=2.»,  or  n  (2) 

1C  \  I  t  IB 

subject  to  internal  stability  of  the  nominal.  Here 

T(P,K)=Pu  *  PuW-PaK)’lPn- 

The  solution  of  this  problem  for  a=2  and  »  is  the  focus  of  Part  2  on  Syn¬ 
thesis  Theory.  The  solution  presented  there  unifies  the  two  approaches  in  a 
common  synthesis  framework.  The  a  s  2  case  was  already  known  and  the 
results  are  simply  a  new  interpretation.  The  o=»  case  had  been  solved  only 
for  special  cases  where  Fn  and  Pzi  are  square.  Also,  the  existing  solutions 
did  not  have  computational  schemes  allowing  their  use  on  even  moderately 
sized  problems.  These  twi  limitations,  especially  the  former,  restricted  the 
applic  .ion  of  the  pioneering  Hm  methods  to  fairly  simple  problems,  such  as 
sensitivity  minimization.  The  new  solution  eliminates  these  two  Limitations. 
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Unfortunately,  this  new  solution  for  the  Hz  and  Hm  suffers  from,  the 
same  limitations  imposed  by  restrictive  assumptions  about  uncertainty  as  do 
the  underlying  analysis  methods.  While  the  S5V  is  a  great  improvement  for 
analysis  (Case  3),  synthesis  for  the  a=fi  case  is  not  yet  fully  solved.  Recalling 
that  (a  may  be  obtained  by  scaling  and  applying  •  a  reasonable  approach 
is  to  'solve" 


(3) 


by  iteratively  solving  for  K  and  D.  With  either  K  or  D  fixed,  the  global 
optimum  in  the  other  variable  may  be  found  using  the  (a  and  Rm  solutions 
described  previously.  Example  designs  have  been  done  and  this  scheme 
seems  to  work  well,  but  global  convergence  is  not  guaranteed.  In  fact,  a 
counterexample  has  been  constructed  where  (3)  reaches  a  local  minimum 
which  is  not  global. 
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2.1.1  Overview  of  Synthesis 

From  the  previous  part  of  these  notes  on  analysis,  we  have  seen  that  the 
synthesis  problem  in  each  case  reduces  to  finding  a  controller  K  which 
achieves  internal  stability  and  solves 

min 

where 

a  P\z 
21  Pa 

and 

Pi{P.K)  *  Pn  +  PM~PaK)'lPzi 

We  will  restrict  our  attention  for  now  to  the'  a=2  and  «  cases  of  (1). 
Recall  that  the  a=(i  case  of  (1)  can  be  converted  to  the  a=«  case  by  scaling. 
The  approach  of  these  notes  is  to  develop  the  o=2  and  -  cases  in  a  parallel 
manner,  emphasizing  their  common  features. 

We  begin  by  considering  the  special  case  of  (1)  where- all  matrices  are 
constants.  This  is  an  interesting  problem  in  its  own  right  and  manages  to 
capture  the  essential  features  of  the  general  problem.  While  the  a=2  case  is 
quite  straightforward,  the  key  step  in  the  solution  of  (1)  for  a  =  «  was  first 
published  in  1982  by  Davis,  Kahan,  and  Weinberger  in  their  important  paper 
on  norm-preserving  dilations. 

The  remainder  of  this  part  of  the  notes  involves  taking  each  step  of  the 
solution  to  (1)  for  the  constant  case  and  generalizing  to  the  case  of  real- 
rational  matrices.  The  difficulty  arises  from  stability /causality  considera¬ 
tions  which  are  not  present  in  the  constant  matrix  case. 


e  Jif  .  pw  e  ^ 


0=2,-,  or  n 


(1) 
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2.1.2  Constant  Matrix  Case 


In  this  section,  we  will  consider  &  special  synthesis  problem  where  ail 
matrices  are  constants.  The  constant  matrix  case  will  allow  us  to  study  the 
synthesis  problem  in  a  simplified  context,  but  one  which  parallels  the 
rational  case. 

For  constant  matrices,  the  norms  reduce  to 


II  P  II.  *  *(P) 

II P  lb  *  (WP)f 

Note  that  these  definitions  are  not  conventional,  but  they  are  convenient  in 
allowing  parallel  development  of  the  constant  and  rational  cases. 

Consider  the  constant  matrix  problem 


min 


Tt(P,K) 


a  =  2," 


(1) 


where 

e  Pi)  e 
and 


P  = 


Pn  P 12 

Pi  l  P  22 


TdP.K)  =  Pu'PiiK(I-PaK)-lPzi' 

Assume  that  P*zP\z  >  0  and  Pz\Pz\  >  0. 

The  first  step  is  to  make  the  substitution  of  variables 

K=Q(I+PzzQ)-\  Q  =  {PSiPi2)-',iQ{Pi:Pii)-'li  (8) 


so 


q  =  (PuPxm-Pziio-'iPuPii  )H 


(3) 


z 


Using  the  linear  fractional  representation  notation, 

K  -  Ti(J,Q) 


(4) 


where 


J\\  J 12 

o  (PStPiz)'* 

?2l  Jtt, 

[P 2\P 21 )  -*  HPzxPzi ) -^Pzl(P^Pu) 

With  this  substitution  we  have 

Ti(P.K)  =  TiiPJi  (/.$)) 

=  Pn  +  [Pu(  Wia)  “»]  JfPzxPZ )  *H^ij 

=  Tu  +  TizQTzi  (5) 


where  the  Ty  are  defined  in  the  obvious  way.  This  parametrization  has  con¬ 
verted  the  nonlinear  problem  in  (1)  to  one  affine  in  the  parameter  Q.  Note 
that  Tufa  =  I  and  Ti\T$i  -  /.  Thus  we  can  find  7j_and  fj^such  that  both 

are  square  and  unitary. 


7*21 

Tu  T'jj  and 

Til 

Since  both  a  -  2  and  ■  norms  are  unitary  invariant 


(T’u  +  T12QT21  a  = 


-1 

Q  0 

J«j 

IL = 

Tu  + 

Ttt  Tjj 

°  o 

Ird 

[r»  t  ,!  t 


Tn 
mrjj 


Q  0 
0  0 


Pn  *  Q  P 12 
P  21  P  22, 


la 


(6) 


where 


m  m  m  mm  mm  m  m m 

T tfi  11  «*  21  *  12*  tlil 

JLTnTZ  T£TuTl 


(7) 
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Thus  the  problem  in  (1)  reduces  to 


Rft  Ilf's* 


i+Q  Fa 

p’ci  ||[  Fzi  R& 


and  soluUon  of  (8)  yields  a  solution  of  (1)  by  solving  (2)  for  K. 
The  a  a  2  case  can  be  solved  immediately  from  (8)  since 


f. 


u *Q  Fu 
F  2i  R22 


\Ru*Q  L  + 


12 


0  R)2 
[Rzi  R& 


Thus 


0) 


O) 


Qa ft  *”•/?!!  *  ~TiaTuT{ i 


and 


Rn+Q  R12 

0  R 12 

Ru  Rzt 

a 

Ru  Rz& 

do) 


The  simplicity  of  the  a=2  case  is  responsible  for  much  of  its  appeal. 

Optimization  in  this  norm  reduces 'to  projection  since  la  is  a  Hilbert  space. 

« 

This  holds  as  well  for  the  rational  matrix' problem. 

The  a=«  case  is.  somewhat  more  complicated  since  Lm  is  not  a  Hilbert 
space  and  the  minimization  in  (8)  cannot  be  solved  by  projection.  For¬ 
tunately.  Lm  arises  as  the  space  of  linear  operators  on  the  Hilbert  space  Lz, 
and  (8)  can  be  treated  as  a  dilation  problem.  The  next  section  focuses  on 
matrix  dilation  problems. 
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2.1.3  Matrix  Dilation.  Problems 


Consider  the  optimization  problem 


min 

x 


X  B 
C  A 


(D 


where  X,B,C,A 

restatement  of 
submatrices  os 


are  constant  matrices  of  compatible  dimensions.  This  is  a 


(2.7)  for  the  a=<*  case.  The  matrix 
indicated  in  the  following  diagram: 


\X  B\ 
C  A 


is  a  dilation  of  its 


c 


X  B 

-4 

B 

C  A 

4- 

A 

d 

d  c 

d 

r 

e 

f  ' 

r 

[C  A 

4- 

\A 

(2) 


In  this  diagram,  c  stands  for  the  operation  of  compression  and  d  stands 
for  dilation.  Compression  is  always  norm  decreasing;  sometimes  dilation 
can  be  made  to  be  norm  preserving.  Norm  preserving  dilations  are  the  focus 
of  this  section. 

The  simplest  matrix  dilation  problem  occurs  when  solving 

(3) 

Although  (3)  is  a  much  simplified  version  of  (l),  we  will  see  that  it  contains  ail 
the  essential  features  of  the  problem.  Letting  ye  denote  the  minimum  norm 


mm 

x 
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in  (3),  it  is  immediate  that  yt  - 1 |A|  |».  The  following  theorem  characterizes  all 
solutions  to  (3). 


Theorem  1:  For  v?  2  yt , 


iff  2  Y  with  Y  1  such  that 


Proof: 


iff 


X*YtfI-AmA)* 


X*X  +  AmA  sc  y2! 
iff 


X'X*(ylI-AmA) 

iff 

||au||s  ||(t *I-AmAp  u||  Vu 
iff 

X  -  Y(yiI-AmA)^  for  some  ||yjp£l 


W 

(5) 


This  theorem  implies  that,  in  general,  (3)  has  more  th«n  one  solution.  This  is 
in  contrast  to  the  a  =  2  case.  The  solution  X  *  0  is  the  central  solution  but 
others  are  possible  unless  A* A  -  7s/.  A  more  restricted  version  of  the 


theorem  is 


3 


Corollary  1:  For  y>yt , 


iff 


X 

A 


*7 


| X&I-A'A)-* 


sS  1. 


(6) 

(7) 


The  corresponding  dual  results  are 


Theorem  1* 


such  that 


For  v?  i  7, 


A  j  |.  <  7  ill 


X  =  tfl-AA*)*Y 


Corollary  1*  For  7>7e 

m 


^tfl-AA^X 


SSI  1 


(8) 

(9) 


(10) 

(ID 


Now,  returning  to  the  problem  in  (1),  let 


y0  =  min 
,a  x 


X  3 
C  A 


(12) 


The  following  theorem  (Parrott)  will  play  e.  central  role  in  the  synthesis 
theory.  The  proof  is  a  straightforward  application  of  Theorem  1  and  1’. 


4 


Proof:  Denote  by  7  the  right  hand  side  of  the  equation  (13).  Clearly,  7*^7 
since  compressions  are  norm  decreasing.  That  ya<7  will  be  shown  by  using 
Theorem  1  and  1*. 


From  Theorem  1  we  have  that  B-Yi^I-A'A)^  for  some  Y  such  that 
Y  «s£l.  Similarly.  Theorem  1*  yields  C-(j*I-AAm)*Z  for  some  Z  with 
Z  |rfl. 

Let  1  a  -YAmZ.  Then 


-YA*Z  Y&I-AmAfl 
[(7 •l-M'fiZ  A  |1- 

[  -4"  (^/-4*4)H| 

1(9*/-Af)»  A  l 


■  ? 

Since 

-A*  (f*/-A“A)»|f  -4  (tV-AA*)*]  tyl  0 

(7 2I-AA*ft  A  II (7*1 -Am A)*  4*  1  =  |o  7 */ 

Thus  */&7t.  so?=7t. 


This  theorem  gives  one  solution  to  (12)  and  an  expression  for  y0.  As  in 
(3),  there  may  be  more  than  one  solution  to  (12),  although  Theorem  2  only 
exhibits  one.  Theorem  3  considers  che  problem  of  parametrizing  ail  solu¬ 
tions.  The  solution  X  -  -YA*Z  is  the  "central"  solution  analogous  to  X  -  0  in 
(3).  The  next  corollary  is  an  alternative  statement  of  Theorem  2  using  the 


form  of  (2.4)  for  the  problem 

7,  =  cn^n  |[fl  +  UQV ||„ 

where  U'U  e  /  and  W  =  / 


5 


(14) 


Corollary  2 


7,  =  max 


IHL.  H 


(15) 


The  following  theorem  (Davie,  Kahan,  and  Weinberger)  parametrizes  all 
solutions  to  (1).  The  proof  ie  omitted,  but  is  similar  to  Theorem  2  and 
involves  application  of  Theorem  1  and  1*. 


Theorem  3  Suppose  7  %  y, .  The  solutions  X  such  that 

If (18> 

are  exactly  those  of  the  form 

X  =  -YA*Z  *  7 {I-YY‘)*W(I-ZaZ)*  (17) 

where  V  is  an  arbitrary  contraction  (||w||«£l)  and  Y  and  Z  solve  the  linear 
equations 

B  =  YtfI-AmA)* 

C  =  (?I-AAm)*Z.  (IS) 
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The  following  corollary  gives  an  alternative  version  of  Theorem  3. 

Corollary  3  For  y  >yt, 

Ilk  s||| 

IImIL*  7 

i iff 

where ■ 

Z  ~  (7*I-AA')-*C 

There  are  many  alternative  characterizations  of  solutions  to  (19), 
although  the  formulas  in  (20)  and  (21)  seem  to  be  the  simplest 

For  the  problem  in  (14),  the  following  equivalences  apply  for  all  7  >  : 


\\R  +  UQv\\m*7 

(22) 

|  RV+UQ  tf*i*]|L*7 

(23) 

||  \R1+VQ  ^]|L*7 

(24) 

—  ot m  _ j  a  —  dt/« 

j  "•  Jfc  r  cuu  a  —  i  *  * 

(y2/-A4°)"^  UQ  ]  jj^JS  1 

(25) 

(  by  Corollary  1* ) 

iff 


(19) 


(20) 


(21) 


? 


Rz+U2Q 


ss  1 


where  Rj  -  (y2/ -AA*)  ~^Ri 


Ui  =  (y3r-AAm)-*U. 


(28) 


To  complete  this,  simply  factor  Uz  to  extract  a  unitary  factor,  and  apply  toe 
dual  of  (22)*{28)  to  (28).  Although  the  formulas  get  messy,  (22)  can  be 
solved  in  this  manner. 

In  each  of  these  cases,  Theorem  2,  Theorem  3,  their  corollaries,  and  the 
solution  described  above,  the  general  case  reduces  almost  immediately  to 
application  of  Theorem  1  or  Corollary  1.  Thus,  when  it  is  convenient,  we  will 
consider  (4)  rather  than  (18)  and  (28)  rather  than  (22).  This  will  simplify  the 
discussion  of  the  rational  case  without  introducing  any  loss  of  generality. 


1 


Z  1.4  Summary  of  Constant  Problem 


The  rational  matrix  problem  in  equation  (1.1)  can  be  solved  in  a  manner 
which  parallels  the  treatment  of  the  constant  case  in  the  last  two  sections. 
This  generalization  is  the  focus  of  the  next  three  chapters  on  synthesis.  To 
reinforce  the  similarity  between  the  constant  and  rational  case,  we  will  now 
review  the  key  steps  from  the  previous  two  sections  and  preview  their  gen¬ 
eralizations  to  the  rational  case. 

Consider  the  diagram  in  Figure  1.  This  summarizes  the  steps  in  the  con¬ 
stant  matrix  problem  (2.1).  The  main  steps  are  as  follows: 


1)  Parametrization:  Make  the  substitution  K-Fi(J,Q)  so  that 

Fi(F.K)  *  Fl(P.Fl(J,Q)) 

*  Fi(T,Q) 

s  Tn  T’sa^T’ai  (l) 

is  affine.  Additionally,  we  want  T^T\z^I  and  TziTzi-I. 


2) 


is  r  i 

Unitary  Invariance:  Find  7^ and  Tj^so  that  \Tl2  Tjj  and 
and  unitary. 


Tz  i 

TjJ 


are  square 


Pre-  ard  post-multiply  by  and 


-]■ 


Ux 


to  yield 


Fu+Q  Bu 

Rz\  Rzz 


where 


(2) 


2 


R  = 

Ris 

Rzi 

Rz 2 

n 2 

_  f 

X 

n. 

M 

(3) 


Recall  that  without  loss  of  generality,  we  may  assume  the  ru?£»0  so 
that  (2)  becomes 

tel 

3)  Projection  /  Dilation;  At  this  point  the  a~ 2  and  o=«  cases  differ.  For 
a- 2,  the  problem  reduces,  by  projection,  to 

nun  \\Rn  *  c||j  (5) 


which  has  the  unique  solution  Q--Ru. 


The  a=«  case  must  be  treated  using  the  matrix  dilation  theory  of  the 
previous  section.  Recall  that,  in  general,  the  solution  is  not  unique. 
From  Theorem  3.1.  all  solutions  to 


for  7  fc  j/?*i  1 1 


(6) 


are  of  the  form 


for  some 
that 


q  =  -  Ytfl-RfM* 


(7) 


«S1.  Corollary  3.1  gave  the  alternative  characterization 


R’l 


7 


for 


IWL 


(9) 


if  and  only  if 


||(*u+$)(72/-*aV?tt)Hi||.  *  L 


(9) 


•It  is  this  latter  characterization  which  will  be  used  in  the  rational  case. 
Note  that  Q--Rn  is  one  solution  to  (8)  and  (8). 

Recovery  of  the  optimal  K:  This  is  obtained  by  simply  computing  K 
from  the  formula  K~F^J,Q)  used  in  step  1)  to  parametrize  the  prob- 


2.1.4  Figure  1 


min  ||W.JO| 


n^n  I] r„  + 


where  F^PJO  =  Pn  +  PM~PaK)"Pti 


TizQTji 


paremetrization 

K~Fi(J,Q) 

where  T^TizsI  failai  s  / 


unitary  invariance 


mp 


Fn—Q  P31 
P\2  F32 


a=m 


a=2 

projection 


dilation 
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2.1.5  Rational  Matrix  Generalization 

The  steps  In  the  rational  ease  closely  parallel  the  constant  case,  as 
shown  in  Figure  1.  Most  of  the  work  in  the  remaining  chapters  is  devoted  to 
generalizing  these  steps  from  constants  to  rationals.  The  source  of  all  the 
difficulty  in  the  rational  case  comes  from  the  requirement  for  internal  stabil¬ 
ity,  or  equivalently,  causality.  Without  this  the  rational  case  would  reduce  to 
the  constant  case  at  each  frequency,  and  could  be  solved  using  the  results  of 
the  previous  two  sections. 

We  will  now  briefly  outline  the  steps  required  to  solve  the  rational  case 
and  preview  the  upcoming  chapters. 

1)  Paramatrization:  Find  J  so  that  the  substitution  K-F^J.Q)  yields 

Fi(P.K)  s  FtiPSdJ.Q)) 

*  Fi(T.Q) 

*  Tu  +  TuQTn  (10) 

with  the  additional  requirement  that  T  e  Hm  and 

Fi{P,K)  internally  stable  (11) 

iff  Q  e  Hm 

This  parametrizes  all  stabilizing  ICa  in  terms  of  a  stable  Q  e  Hm  in  addi¬ 
tion  to  providing  an  affine  parametrization  of  all  stable  Fi(P,K).  This 
parametrization  (Youla)  is  developed  in  Chapter  2  on  Stabilization. 

A  further  requirement  is  that  T&  and  Tzi  be  inner,  that  is  TxtTiz-^  and 
Tz\Tzi  -I.  Methods  for  obtaining  the  particular  parametnzations  which 
achieve  this  are  developed  in  Chapter  3  on  Factorization. 
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1 

Tn 

T 12  7j|  and 

Mm 

[TjJ 

are 


2)  Unitary  Invariance:  Find  T and  7j_  so  that 

squart  and  inner  (also  Chapter  3).  Then  pre-  and  post-multiply  by 


f  u 

Tt  i 

|r12  fjj  and 

to  yield 


f- 


n+9  R ia 

R  ai  R72 


(12) 


whore 


R 


Again,  to  simplify  the  presentation  suppose  that  so  that  (12) 

becomes 


3) 


Projection  /  Dilation:  At  this  point  the  a=2  and  a=» 
For  a=2.  the  problem  reduces,  by  projection,  to 


fSg.  li*»  *  «IL 


(13) 


cases  again  differ. 


(14) 


But  since  R\\Z  Lm  ,  Q=-R\\  would  not  correspond  to  a  stable  solution. 
The  unique  solution  is  yet  another  projection 

QofA  -  Pnt(Rn)  (15) 

where  Pnx  denotes  projection  onto  Hz-  When  viewed  appropriately,  these 
two  projections  can  be  seen  as  a  single  projection  onto  a  subspace  of 
Lz(j  R  t'™1). 


The  o=“  case  is  again  treated  as  a  dilation  problem.  Since,  genericaliy 


3 


min 

<  Hm 


*u+Q 

*n 


>  7  I 


(18) 


it  is  convenient  to  use  the  characterization  in  Corollary  3.1  and  (8)-{9). 
Recall  that 


iff 


Rn+Q 

*21 


*7 


for  7  >7 


(IV) 


li(/?u+<?)(7a/-^2i)"H|L 56 1  '  (is) 


The  key  to  proceeding  in  the  rational  case  is  to  find  M  G  RHm  such  that 
U~l  G  RHm  and  =(7*1  -RziRu)-  If  we  use  the  symbol  {t3! -*21*21)^ 

to  denote  this  M,  then  (18)  makes  sense  in  the  rational  case.  Finding  M 
involves  spectral  factorization  and  is  treated  in  Chapter  3. 


Given  U  6  RHm  with  the  desired  properties,  (18)  reduces  to 

\g  +  9||-  *  1  (19) 

where  C~RnH~l  G  RHm  and  §= Wl.  Solving  (19)  for  Q  G  RHm  solves 
(18)  for  Q  G  RH^  Note  that  is  in  RHm  if  §  is,  since  M  e  RHm  by 

construction. 


The  final  step  in  the  rational  case  then  involves  solving  (19)  for  §  G  RHm- 
This  is  a  standard  mathematical  problem  of  approximating  an  Lm  matrix 
by  an  Hm  matrix  This  turns  out  to  be  yet  another  dilation  problem  but 
in  a  somewhat  different  guise  than  those  treated  in  the  constant  case. 
The  solution  of  (19)  is  the  focus  of  Chapter  4  on  Best  Approximation. 

4)  Recovery  of  the  optimal  K:  Just  as  in  the  constant  case  ~Fi(J,  Q„pt ). 
This  Kopt  will  stabilize  Fi(P,K ^,t)  since  the  parametnzation  in  Step  i) 
insured  that  Q  stable  lead  to  internal  stabliiity  of 


2.1.5  Figure  1 


min 


|  (P, K)  |  |a  internal  stability 


parametmation 

K=Ft(J.Q) 


0**$  Tu  +  ^12^*21  |a  "fcere  T*2T\2~I  TziTzi-I 


f 


unitary  invariance 


A 


|i?n“9  Rai 
R 12  R  22 
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2.2.2  Internal  Stability 

In  this  section  P  and  K  are  fixed  proper  transfer  matrices.  The  block 
diagram  of  Figure  1  represents  the  two  equations 


Partition  P  accordingly; 


a-4- 

Pz\  P&  1 


It  is  convenient  to  introduce  two  fictitious  external  signals,  and  tua.  as  in 
Figure  la. 

Suppose  the  signals  v  .u,.  and  Wz  are  specified  and  that  u  in  Figure  la  is 
well-defined.  Then  so  are  e  and  y.  Thus  it  makes  sense  to  define  the  system 

v 

in  Figure  la  to  be  wall  -posed  provided  the  transfer  matrix  from  to  u 

w2 

exists  and  is  proper. 

Lemma  1.  The  system  is  well-posed  if  and  only  if 

/  -  KWP&i")  is  invertible.  (3) 


Proof.  Figure  la  implies  the  equations 
u  ~  wx  -  KY  *  Kvjz 

y  s  Pz\v  +  Pziu 

and  these  in  turn  imply  that 

(I-KP)u  -  w  i  *  KP ziV  +  Kwz. 

Thus  well-posedness  if  equivalent  to  the  condition  that  ( I-KP)~l  exists  and  is 


proper. 


2 


QED 

It  is  straightforward  to  show  that  (3)  ia  equivalent  to  either  of  the  follow¬ 
ing  two  conditions: 


/  -*(“) 
—Pa(m)  I 


is  invertible : 


I  -  Pa(m)K(m)  is  invertible  . 


(4) 

(5) 


The  well-posedness  condition  is  simple  to  state  in  terms  of  state-space 
realizations.  Introduce  minimal  realizations  of  P  and  K: 


P  = 


A  B\  B% 
Ci  Bn  Diz 
Ca  D3l  Dzi 


K 


A  § 

t  tr 


(0) 

(?) 


Note  that  the  partition  in  (6)  corresponds  to  that  in  (2),  Le., 


Pa 


A  B) 


(8) 


Then  Pzzi^-Dza  and  K(m)-S,  so  for  example,  from  (4)  well-posedness  is 
equivalent  to  the  condition  that 


is  invertible. 


0) 


Well-posedness  will  be  assumed  for  the  rest  of  this  section.  Let  x  and  S 
denote  the  state  vectors  for  P  and  K  respectively,  and  write  the  system 
oCjUntions  in  5*i**urs  i  v  to  2oro  2nd.  _  i^no^sd! 


2  -  Ax  *  Bau 
y  -  Cjx  -  D&u- 


(10a) 

(10b) 
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AS  +  By 

(10c) 

Bs  +  By. 

(lOd) 

The  system  of  Figure  1  Is  internally  stable  provided  the  origin  (x,S)  =  (0,0) 
is  asymptotically  stable.  To  get  a  concrete  characterization  of  internal  sta¬ 
bility.  solve  equations  (10b)  and  (lOd)  foru  andy: 


u  I  -B 

-i 

o  ?][* 

y  *  [-£»  / 

c2  o][s 

(Note  that  the  inverse  exists  from  (9)).  Now  substitute  this  into  (10a)  and 
(10c)  to  get 


jeL* 

dt  £ 


where 


A  0 

Bi  0 

I  -B 

-i 

0  t 

P  2. 

o 

-Dzi  I 

C3  0 

Thus  internal  stability  is  equivalent  to  the  condition  that  A  has  all  its  eigen¬ 
values  in  the  open  left  half-plane. 

It  is  routine  to  verify  that  the  above  definition  of  internal  stability 
depends  only  on  P  and  K ,  not  specific  realizations  of  them.  The  following 
result  is  standard. 

Lemma  2.  Consider  a  minimal  realization  of  P  as  in  (6).  There  exists  a 
proper  K  achieving  internal  stability  if!  (A,32)  is  stabilizable  and  {C2A)  is 
detectable. 

The  latter  stabilizability  and  detectability  conditions  are  assumed 
throughout  the  remainder  of  this  chapter.  Since 


4 


P-a 


A 

Bz 

Ct 

Dzz 

(U) 


equations  (10)  constitute  a  state-space  representation  of  the  system  in  Fig¬ 
ure  Z  Although  the  realization  in  (11)  is  not  necessarily  minimal.  It  is  stabil- 
izable  and  detectable,  and  these  are  enough  to  yield  the  following  result 

Lemma  3.  The  system  in  Figure  1  is  internally  stable  if!  the  one  in  Figure  2 
is. 


The  next  section  contains  a  parametrization  of  all  K" s  which  achieve 
internal  stability  for  the  system  in  Figure  Z  To  simplify  notation,  define 

G  Pzz  1  B  :s  Bz  •  C Cz ,  D  :*  D33. . 

Then  ( A,B )  is  stabilizable,  (C.A)  is  detectable,  and  the  system  under  study  is 
that  in  Figure  3. 


The  above  notion  of  internal  stability  is  defined  in  terms  of  state-space 
realizations  of  G  and  K.  It  is  important  and  useful  to  characterize  internal 
stability  from  an  input/output  point  of  view.  For  this,  consider  the  feedback 
system  in  Figure  4.  This  system  is  described  by. 


Now  it  is  intuitively  clear  th&i.  if  the  system  in  Figure  4  is  internally  stable 
then  for  all  bounded  inputs  (t^.i/j),  the  outputs  (ct,^)  are  also  bounded.  The 
following  lemma  shows  that  this  idea  lends  to  an  input/output  characteriza- 
Mrtie  stability. 

Lemma  4.  The  system  in  Figure  4  is  intsrnaiiy  stable  if  and  only  if  {I-GK) 
is  invertible  and  the  transfer  matrix 
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I  -K 

-l 

r+K(I~GlTlG  K(I-GK)-1 

-G  I 

7 

1 

• 

(ID 


between  (vj.Vj)  and  (*i.*a)  belongs  to  RH~ 


Proof.  As  above  let  ( A,B,C,D )  and  (A>B,B,3)  be  stabilizable  and  detect¬ 
able  realizations  of  G  and  K  respectively.  Then  the  si'  te-apace  equations  for 
the  system  in  Figure  4  are: 


\B  C  «ij 
0  B  Bgl 


D 

0 


0 

3 


•  i 

8  a. 


St  =  v:  +  yj.  »t  =  v-i  +  y  i. 


The  last  two  equations  can  be  rewritten  as 


/  ~d\ 

[o  t 

*1 

«i 

rD  /] 

8* 

"  [C  0 

H 

+ 

Now  suppose  that  this  system  is  internally  stable.  Then  (7)  implies  that 
-  {I”GK)(.<*)  is  invertible.  Hence  (I-GK)  is  invertible.  Further, 
since  the  eigenvalues  of 


U  o. 

B  o] 

I  -D 

-l 

0  C 

Ml  * 

o  .9; 

-D  I  _ 

C  0 

are  in  the  open  left  half  plane,  it  follows  that  the  transfer  matrix  in  (10)  from 
(Vi.v-B  to  (8  i.Ug)  is  in  RH*. 

Conversely,  supppose  that  ( I-GK)  is  invertible  and  the  trrnsfer  matrix 
in  (i0)  is  in  Rfla.  Then,  in  particular,  (/— Ga)~1  is  prcpei  which  implies  that 
(I—GK)(eB)  =  {I-DB)  is  invertible.  Therefore 


/  -3 

-2?  / 


6 


is  nonsingular.  Now  routine  transfer  function  calculations  give, 


/  -j 3  *i 

0  t 

B  oil 

-D  I  jej]  = 

I  + 

C  0 

(S/-A)"1 

0 

Since  the  transfer  matrix  from  (t'j.vj)  to  («i.*8)  belongs  to  RH*.,  it  follows 
that 


0  3 

B  0 

C  0 

c *i-a r1 

o  B 

belongs  to  RHm.  Finally,  since  ( A.B.C )  and  (2,3.3)  are  stabiiizable  and 
detectable,  it  follows  that  the  eigenvalues  of  A  are  in  the  open  left  half  plane. 


QED 

We  note  that  to  check  internal  stability  it  is  necessary  (r.nd  sufficient)  to 
check  that  each  of  the  four  transfer  matrices  in  (J.1)  ere  in  RHm.  It  is  not 
difficult  to  construct  examples  of  C?  and  K  such  that  any  three  of  the  four 
transfer  matrices  in  ('1  i)  arc  in  RHm  while  the  fourth  one  is  unstable. 
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2.2.3  Parametrization  of  All  Stabilizing  Controllers 


Two  matrices  N,£i  e  RHm  with  the  same  number .  of  columns  are 


Ur 


right  -coprima  if  the  combined  matrix  has  a  left  inverse  in  RH *  That  is 


there  exists  X,Y  Z  RHm  such  that  XM  +  YN  - 1.  This  is  often  called  a 
Bezout  or  Diophantine  equation.  An  alternative  definition  is  that  two 
matrices  in  RHm  are  right-coprime  if  every  common  right  divisor  in  RHm  Is 
invertible  in  RHm-  This  can  be  shown  to  be  equivalent,  to  the  above  definition 
in  terms  of  a  left  inverse,  but  we  will  not  use  this  fact. 


It  is  a  fact  that  every  G  Z  (proper,  real-rational)  has  a  right-coprime 
factorization  G  -  NH~l  where  N,H  C  RHm  are  right  coprime.  Similarly, 
there  exist  left  coprime  factorizations  (lcf),  defined  in  the  obvious  way  by 
duality.  The  proof  of  the  existence  of  such  coprime  factorizations  will  be 
given  in  the  next  section  with  explicit  realizations  for  the  factorizations.  In 
this  section,  we  will  see  how  these  factorizations  can  be  used  to  obtain  a 
parametrization  of  all  stabilizing  controllers. 

Begin  with  refs  and  lcf  s  of  G  and  K  in  Figure  4: 

G  =  NM~X  =  M~XN  (1) 

K  -  irrx  s  r‘xu  (2) 


Lemma  1.  Consider  the  system  in  Figure  4.  The  following  conditions  are 
equivalent: 


1. 

Z 


The  feedback  system  is  internally  stable. 

\u  rjl 


\ff  v\ 


is  Invertible  in  RHm. 


3. 


I 

V  - U 
\-N  M 


is  invertible  in  RH, 


2 


Proof:  As  we  saw  in  Lamma  2.3  of  the  last  section,  internal  stability  is 
equivalent  to  the  condition  that 


/ 

-ff 


€  RE. 


or.  equivalently, 


i  I  X 
iff  / 


-i 


e  RHm. 


(3) 


Now 


/  K 

/ 

UVl 

G  I 

NM~l 

/ 

. 

U  V 

kl 

o  ! 

“ 

N  V 

[  o  r-lj 

so  that 


I  K 

'*  [wo 

u  u 

G  I 

*  [o  V 

N  V 

Since  the  matrices 


u  o]  \m  d\ 

0  ij*  [N  Vj 


are  right*coprime,  (3)  holds  iff 


Uf  U 
LN  V 


-i 


e  RE*. 


This  proves  the  equivalence  of  conditions  1  and  2.  The  equivalence  of  1  and  3 
is  proved  similarly. 


3 


We  shall  see  in  the  next  sectjon  how  to  And  explicit  realizations  for 

N,  M,  N,  M,  Ut,  Vt,  Ut,  V,  and.such  that  (l)  holds  and 

1 

V ,  -ut 
[-N  U 

The  above  lemma  says  that 

a;  S  u,  ir1  = 


]m  V, 

I  0 

U  V, 

0  / 

(5) 


then  qualifies  a  particular  controller  achieving  internal  stability.  All  stabiliz¬ 
ing  controllers  can  be  expressed  in  terms  of  K,  and  a  parameter  Q,  as  shown 
in  the  following: 

Theorem  1.  The  set  of  all  proper  controllers  achieving  internal  stability  is 
parametrized  by  the  formula 

K*K,  +  V^Q(I+V?'NQ)-lKx  (8) 

where  Q  ranges  over  RHm  such  that  (/  +  VTlNQ)(-)  Is  invertible. 

Proof:  Assume  K  has  the  form  indicated. 


Define 


U  ft  U0  *  MQ,  V  £  V9  +  NQ 
U  £  Ua  +  QM,  Vi  V0  +  QN 


then 


V 

-V 

u  u 

V0  +QN  -(Ug+QM) 

-N 

M 

N  V 

-N  M 

M  U, 

N  V0+NQ\ 


\i  -Qp,  -Sr. }\a  %][/  d 
“  [0  /  II-/7  M  j[AT  Vg  |[0  I\ 


I  -Q 

I  Q 

0  / 

o  / 

from  (5) 
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I  0 

=  1°  1 

Thus  K  achieves  internal  stability  by  lemma  1. 


(7) 


Conversely,  suppose  K  is  proper  and  it  achieves  internal  stability.  Intro¬ 
duce  ref  and  lef  of  K  as  in  (2). 


Then  by  lemma  1,  Z  k  HV  -  NU  \a  invertible  in  RH*.  Define  Q  by  the 
equation 


U,  +  UQ*  U2T\ 


(8) 


Q  =  -  Ut) 

Then 

V,  +  NQ  s  V,  +  NH-\UZ~X  -  Ut) 
a  V,  +  M~lB{UZ~x  -  Ut) 

*  S~l(Svt  -  NU„  +  NUZ~X) 
=  M~X(I  +  jVC/Z"1) 

*  ir'C?  + 

=  M~XMVZ~X 

*  VZ'1 


Thus, 


/r  =  t/y‘ 


-  (U'+MQKV,+NQ)-x 


—  rf  tr 

-  Wp 


/  f/ 


rr  t/—  1  a?\  n/  r 
wp  ry  j* 


t/-1  *rn\— 1 
ro  '• 


O) 


from  (1) 

from  (5) 


(10) 


from  (prelim?  ).  Then,  since 

{M  -  V,  V?XN)  =  (U  -  V~'U0  V)  =  V7l(K0#  -  C7„JV)  =  V~l 


/i  i  \ 


5 


we  have  that 

K  s  U'V'  +  V,lQ(I  +  V?lNQ)-lV?1.  (12) 

To  see  that  Q  belongs  to  RHm,  observe  first  from  (9)  and  (10)  that  UQ 
and  NQ  both  do.  Right-coprim eness  of  N  and  U  then  implies  that  Q  e  RH*. 

Finally,  since  V  and  Z  evaluated  at  s==*  are  both  invertible,  so  is 
Va  +  NQ,  from  (10),  hence  so  is  I  +  VfxNQ. 

QED 


Define  the  rational  matrix 


J  £ 


K.  V,~l 
VT1  -V?lN 


(13) 


and  consider  a  controller  K  given  by  formula  (S).  Then  the  controller  equa¬ 
tion 


u 


-  Ti(J.Q) y 

-  \k,  +  V-'Q(I  +  V,-'NQ)‘lV-' 


V 


is  equivalent  to  the  triple  of  equations 

u  =  K,y  t 

*i  =  V-h,  -  YrlNy, 

Vi  ~  Qut 

The  biock  diagram  corresponding  tc  this  triple  is  in  Figure  5.  We  con¬ 
clude  that  every  stabilizing  controller  can  be  represented  as  K  -  ¥i(J,Q),  as 
in  Figure  5,  for  some  parameter  Q,  which  is  constrained  only  to  be  stable  and 
proper  and  to  make  K  proper. 

The  next  section  gives  an  explicit  state-space  realization  of  one  choice  of 


the  interconnection  matrix  J. 


1 


2.2.4  Realization  of  J  : 


Recall  from  Section  2  that  we  have 


G* 


A  B 
CD  • 


where  ( A,B )  is  atabilizable  and  ( C,A )  is  detectable.  To  obtain  a  right- 
coprime  factorization  of  G,  choose  a  matrix  F  such  that  A+BF  is  stable. 
Lamma  1.  A  stabilizing  state  feedback  F  yields  ref  G  =  N£i~l  where 


M 

N 


A+BF 

B 

:= 

F 

I 

C+DF 

D 

Proof:  That  G  -  NM~X  follows  from: 


A 

B 

A+BF 

B 

C 

D 

F 

Pi 

A  BF 

B 

s 

0  A+BF 

B 

[C  DF 

D 

A+BF  BF 

B 

sr 

0  A 

0 

C+DF  DF 

D 

A+BF  5 
C+DF  D 


(cascade  of  two  systems) 


(by  change  of  basis  in  the  state-space) 


(deletion  of  uncontrollable  part) 


(1) 


-N. 

That  N  and  if  are  right-coprime  will  follow  from  (3)  below. 


QED 
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A+BF 

B 

:= 

F 

I 

C+DF 

D. 

(la) 


is  also  a  realization  of  an  ref  of  G. 

By  duality,  to  get  a  left-coprime  factorization  of  G,  take  H  such  that 
A+HC  is  stable. 


Lemma  1*.  A  stabilizing  output  injection  H  yields  lef  G-M^N  where 


A+HC\ 

H  B+HD 

C 

I  D 

(2) 


The  next  step  is  to  specify  U0 ,  V0 ,  U9t  Vt  to  satisfy 


'  »*»  ' 
v.  -ua 

M 

u. 

I  0 

~N  M 

N 

V, 

0  I 

The  idea  behind  the  choice  of  these  matrices  is  as  follows.  Using 
observer  theory,  find  a  controller  K,  achieving  internal  stability.  Perform 
factorizations 

K>  =  U0Ve-i  =  V0-'U0 

analogous  to  the  onas  just  performed  on  G.  Then  Lemma  3.1  implies  that  the 
left-hand  side  of  (3)  must  be  invertible  in  RH We  shall  see  that,  in  fact,  (3) 
is  satisfied. 

The  equations  for  K,  are 

£  =  A2  +  Bu  +  H(C£  -  £u_-  y ) 
u  s/2. 


that  is. 


3 


* 


Define 

A  :=  A  +  BF  +  HC  *  HDF  .  5  s« 
£:«/•  ,  J5:*0 

f:*C  +  DF,R:*-{B  +  HD). 
Following  (l)  and  (2),  define 


Using  the  above  definitions  we  have  that 


end  the  following  theorem  holds. 


Theorem  2. 

Equation  (3)  is  satisfied. 

Proof:  Verification  of  (3)  is  immediate  using  (7),  (8).  and  the  inversion  for¬ 
mula  for  systems  (prelim). 

A  realization  of  J  is  now  immediate.  Substitution  of  (1).  (4).  (5), 
and  (6)  into  (3.13)  leads  after  simplification  to 


4 


A+BF+HC+HDF 

-H  B+HD 

J  = 

F 

C  / 

-(C+DF) 

I  -D  , 

Let's  recap.  We  began  with  tbs  (stabilizable  and  detectable)  realization 


A 

B 

C 

D 

We  chose  F  and  H  so  that  A+BF  and  A+HC  were  stable.  Define  J  by  (9). 
Then  the  proper  JTs  achieving  internal  stability  are  precisely  those 
representable  as  in  Figure  5.  where  Q  €  RHm  and 

/  +  DQi-f  is  invertible 
(The  last  condition  is  equivalent  to  the  one 

(I  *  VTlNQ)(m)  is  invertible 
which  is  required  as  per  Theorem  1). 

This  representation  result  has  an  interesting  interpretation  :  every 
internal  stabilization  amounts  to  adding  stable  dynamics  to  the  plant  and 
then  stabilizing  the  extended  plant  by  means  of  an  observer.  The  precise 
statement  is  as  follows;  for  simplicity  of  the  formulas,  only  the  case  of 
strictly  proper  G  and  K  is  treated. 

Theorem  2. 


Assume  G  and  K  are  strictly  proper  and  the  system  in  Figure  3  is  inter¬ 
nally  stable.  Then  G  can  be  embedded  in  a  system 


^  | 

B . 

c, \ 

0 

where 


(10) 


5 


and  Aa  is  stable,  such  that  K  has  the  form 


where  A,  +  B,F,  and  At  +  H,  C,  are  stable. 

Proof.  K  is  representable  as  in  Figure  5  for  some  Q  in  RHm.  For  K  to  be 
strictly  proper,  so  must  Q  be  (see  (3.8)).  Talce  a  minimal  realization  of  Q: 

Q  s  [q,  0  • 

Since  Q  e  RHm,  A*  is  stable.  Let  x  and  s,  denote  state  vectors  for  /  and  Q 
respectively,  and  write  the  equations  for  the  system  in  Figure  5  : 

x  =  (^4  +  BF  HC)x  -By  +  By  j 
u  =  Fx  *  Vi 
M  =  -Cx  *  y 
z„  =  AqX9  +  Buu: 

V\  =  c;za 

These  equations  yield 

*«  =  (A  *  B,F.  +  H,Ct)xt  -  H,y 

u  -  F,x,  , 
where 


»  •  «  - _ _ t  _  /4 

ana  At ,  nc ,  *~i  arc  aa  lu  (iu; , 


QED 
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2.2.5  Closed-Loop  Transfer  Matrix 

Theorem  1  provides  a  parametrization,  in  terms  of  Q.  of  all  proper  IT s 
which  achieve  internal  stability  in  Flours  1.  The  goal  in  thia  section  is  to 
express  the  transfer  matrix  from  v  to  s  in  terms  of  Q. 


A  stabilizing  K  is  representable  as  in  Figure  2.  Substitution  of  the  block 
diagram  in  Figure  5  into  that  in  Figure  i  leads  to  the  one  in  Figure  8.  Elimi¬ 
nation  of  the  signals  u  and  y  leads  to  Figure  7  £er  a  suitable  transfer  matrix 
T.  Thus  all  closed-loop  transfer  matrices  are  rspre?ectabk  m  in  Figure  7.  It 
remains  to  give  a  realization  of  T. 

ft 

We  must  first  put  back  the  original  notation  which,  was  simplified  at  tbs 
end  of  Section  2.  Let 


A 

8 x  B 

p  = 

c , 

Dxx  8 is 

Dz\  Bzl 

be  a  minimal  realization  of  P,  and  choose  F  and  B  so  that  A+BjF  and 
A+HC%  are  etable. 


Lemma.  4. 


A+BzF 

-BZF 

Bx 

Bs 

0 

MEC-i 

Bx+HDzt 

q_ 

"‘x-DnF 

-DxzF 

Du 

&12 

0 

cz 

Dzi 

0 

Proof.  With  the  original  notation  we  have  from  (..4.9)  that 


U+BtF+ffCt+HDisF 

“~H  Bi+HD%2 

J  » 

F 

0  / 

-{c^DkF) 

I  -Dzz 

(2) 


(3) 


Partition  J  and  T  accordingly: 


3 


A+S9F 

Bz 

Biz 

\A+ffCz\ 

Bx^HB  8i] 

m 

Am  \ 

Tzt  K  0. 


la  Figure  7  the  governing  equations  are  therefore 

*  a  Tnv  +  T\zVx 
Uj  =  Tziv 
Vi  51 


(5b) 

(5c) 


so  that 

*  *  (  *i*  7 is QT ti)v- 

In  summary,  we  have 

Tkscmxi  3,  The  set  o i  all  closed-loop  transfer  matrices  from  v  to  e  achiev¬ 
able  by  an- internally  stabilising  proper  controller  is  equal  to 

ru  *  TwQTto  :  Q  £  RHm  1  •>  VztQ(m)  invertible  . 

Tbs  important  points  to  note  are  that  the  closed-loop  transfer  matrix  is 
simply  on  aflina  function  of  the  controller  parameter  matrix  Q  end  that  the 
coefficient  matrices  have  vary  simple  realizations,  namely,  as  in  (5). 


» 


3 


Ta 

Tn 

Tzt 


[  a+b2f 

Bz 

fii  +  DitF 

A+HCt 

B\  +  HDz\ 

Ct 

Dz\ 

0. 


In  Figure  7  the  governing  equations  are  therefore 

*  s  Tnv  +  T12Vi 
uj  =  f21v 
Vi  *  Q* i. 
so  that 

t  *  (Tu  +  TuQTuJv. 


(5b) 

(5c) 


In  summary,  we  have 

Theorem  3.  The  set  of  all  closed-loop  transfer  matrices  from  v  to  «  achiev¬ 
able  by  an- internally  stabilizing  proper  controller  is  equal  to 

Tu  +  TviQTzi  5  Q  e  /  +  *•*$(")  invertible  . 

The  important  points  to  note  are  that  the  closed-loop  transfer  matrix  is 
simply  an  afline  function  of  the  controller  parameter  matrix  Q  and  that  the 
coefficient  matrices  7y  have  vary  simple  realizations,  namely,  as  in  (5). 


1 


2.3.2  Riccati  Equations  and  Factorizations 


Consider  the  Algebraic  Riccati  Equation, 

'  E'X  +  XE  -  XWX  *  Q-0  (ARE) 

where 

E.W.Q  e  K"*",  W  =  W'  i  0  and  9  =  9' 
with  the  associatrd  Hamiltonian  matrix 


E  -V 
-Q  -E' 


(Hamiltonian) 


The  following  theorem  and  corollary  characterizes  the  relationship 
between  spectral  factorization,  Riccati  equations,  and  decomposition  of  Ham¬ 
iltonians. 


Theorem  1  Let  A,B.P,S,R  be  matrices  of  compatible  dimensions  such  that 
P-P' ,  R-R’> 0.  with  (A.5)  stabilizable  and  (P^A)  detectable.  Then  the  fol¬ 
lowing  statements  are  equivalent. 


a)  The  parahermitian  rational  matrix 

r(s)  =  \b'(-si-at1  /][$.  *J 


rl ]p  s 

(sJ-A)~lB 

]  S'  R 

I 

satisfies 


T(j  o)  >  0  for  all  0  s:  •• 


b)  For  E-A-BR~lS' ,  W-BR~lB'  and  Q=P-SR~lS\  there  exists  a  unique 
real  X-X‘  such  that 


2 


E'X  +  XE  -  ^9  =  0 

and  E-BR~xB'X  is  stable, 
c)  The  Kamn \.jnian  matrix 

A-BR~lS'  -BR~lB' 

**  s  I-p+sj?-^  -iA-BR-'sy, 

has  no  j  a*axis  eigenvalues. 

Corollary  1  If  the  conditions  in  Thee*,  am  1  are  satisfied  then  BM  €  Rp 
such  that  M~l  c  RHm  and 

T  =  H'RM 

A  particular  realization  of  one  such  U  is  , 

li  s 

where  F~-R~X{S'* B'X). 

Proof:  (a)-*(c)  Let 


Then  Ah  ~  A  -b5~1C.  Suppose  Ah  has  an  eigenvalue  on-the  j'a-axis.  Then 
3a# ,  x#=(xj'1s?’)‘  such  that 

AffX,  *  ;«#x,  (3) 

or  - 


(jueI’-A)xl  =  -BR~l(S'Xi^  B'zz) 
(jy#/+/.')xj  =  —(P—SR~xS')x\’*SR~xB'Zz 


(4) 

(5) 


3 


Suppose 

0  =  Cz,  =  S'Zi  +  B’zs 

(6) 

then  from  (4)  and  (5) 

(ju*J  “  ■A)*i  *  o  ■ 

(?) 

(jOtI  +  A‘)Z2  -  -PZi 

(8) 

Since  (7)  implies  z*(jutI+A')=Q,  from  (S)  we  have  This  implies, 

along  with  (7)  that  ( P,A )  is  not  detectable.  Hence  Cx,# 0.  Now  Lemma 
0.2.3.8  implies  that  there  exists  u,  w0  such  that  T(j  u)u,  =0.  This  contradicts 
the  hypothesis  that  T{ju)>Q.  Hence  (a)  «♦  (e). 

(b) -*  (a)  Suppose  sX=X  such  that  E-BR~lBX  -  A—BR~l(S'  +  B’X)  is 
stable.  Let  F--R~l{S'*B'X)  and 


A 

B 

-F 

I 

It  is  easily  verified  by  use  of  tlje  Riccati  equation  for  X  and  routine  algebra 
that  f  *  U*RU  so 


rl<s)  =  U-\s)R-\U'{s)Yx 


Now 


= 


A+BF 

B 

F 

7 

So  M~l  e  RHm.  Thus  r-1  e  RLm  and  for  all  0 


.  r-Uif. >\  =  u-iNr.Apnrf-iuW-i  s  n 

«  —  /  —■  v  ”/*’  \*~  %  0  “//  *■  • 


Hence  T(ju)  >  0  and  (b)  -*  (a). 


(c)  -»  (b)  -This  is  proven  as  part  of  Theorem  3.1  in  the  next  section. 
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QED 


The  next  section  focuses  on  the  solution  of  the  Riccati  equation  end  com* 
pletes  the  proof  of  Theorem  1. 


i 


2.3.3  Solution  of  the  Algebraic  Riceati  Equation : 

Consider  once  again  the  Algebraic  Riccati  Equation, 

E7X  *  XE  -  XWX  +  Q  =  0  (ARE) 

•where 

e.W.Q  e  IT*".  W  S  WT  *  0  and  Q  =  C?r 
•with  the  associated  Hamiltonian  matrix 


E  -W 

-<?  -£r 


(Hamiltonian) 


Our  main  interest  is  to  find  the  unique  real  symmetric  stabilizing  solution 
such  that  the  matrix  (E  —  VX)  is  asymptotically  stable.  For  simplicity  we 
•will  use  "solution"  of  the  ARE  to  mean  a  real  symmetric  one.  The  ARE  con¬ 
sidered  here  is  more  general  than  the  ARE  which  arises  in  linear  quadratic 
optimal  control  and  Kalman-Bucy  Altering  theory  in  that  there  is  no  assump¬ 
tion  on  the  deAniteness  of  the  matrix  Q. 


An  important  property  of  the  Hamiltonian  matrix  As  is  that  the  distribu¬ 
tion  of  its  eigenvalues  (denoted  as  A(A# ))  is  symmetric  with  respect  to  both 
the  real  and  imaginary  axes,  i.e ,  if  X  €  MAh)  with  multiplicity  A;,  so  is 
X,  -X,  and  -X.  Therefore,  A  can  be  partitioned  as  Aj  and  Ag  so  that 
X  e  Aj  with  multiplicity  k  implies  that  X  e  Aj  and  -X  ,-X  e  Ag  all  with 
the  same  multiplicity. 


One  connection  between  the  ARE  and  As  can  be  seen  by  assuming  that  X 


is  a  solution  to  ARE  and  conjugating  As  by 


/  0 
\-X  I 


to  yield 


z 


I  0 
-X  I 


E  -W  \l  0 

-q  - et\  [ar  / 


E-WX  —W 

-{. ErX  +  XE-XWX+Q )  -(. E-WX)T\ 


-wx  - w 

0  -(E-WX)T 


This  puts  Aff  in  block  upper  triangular  form  and  clearly  exhibits  a  particular 
partitioning  of  the  eigenvalues  of  A#  with  respect  to  the  Imaginary  axis.  For 
example,  if  E-WX  has  all  its  eigenvalues  in  C_,  then  —(E—WX)7  has  all  its 
poles  in  C+.  Thus,  the  solution  of  /iKE  which  stabilizes  E-WX  yi  lid  a  decom¬ 
position  of  Aj{  into  stable  and  unstable  parts.* 


This  section  will  explore  the  coudiuons  under  which  the  desired  solution 
of  ARE  exists.  There  is  a  considerable  literature  addressing  the  theory  of 
ARE.  and  it  is  not  the  purpose  of  this  report  to  give  a  detailed  treatment  of 
this  subject  We  will  simply  present  and  prove  the  results  which  are  relevant 
to  the  factorization  theorems  in  this  report 


Now,  we  are  going  to  state  the  main  theorem  of  this  section  which  gives 
the  necessary  and  sufficient  conditions  for  the  existence  of  a  unique  stabiliz¬ 
ing  solution  of  (ARE).  Without  loss  of  generality,  we  will  assume  that 
W  =  GGt, 


Theorem  1 : 

The  stabilizability  of  ( E,C )  and  Re[\i(%)3  *  0  ( V  is  1,2, . 2 n)  is 

necessary  as  well  as  sufficient  for  the  existence  of  a  unique  stabilizing  solu¬ 
tion  of  (ARE). 
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Remark : 

■  The  unique  stabilizing  solution  of  Theorem  1  will  be  denoted  by  Ric  (4#). 
Note  that  this  theorem  is  more  general  than  Theorem  2.1  from  the  previous 
section  since  no  detectability  assumptions  are  made.  The  following  theorem 
will  play  un  important  role  in  the  next  section  in  obtaining  complementary 
inner  factors. 

Theorem  2: 

If  Q  -  HrH  i  0  in  (ARE)  and  X  is  its  solution,  then  Ker(X)  c  YJsc{H). 


The  remainder  of  this  section  is  devoted  to  proofs  of  Theorems  1  and  2. 
Lemma  1 :  ( Potter.  Hartensson  ) 


Let  the  columns  of  the  matrix 


€  JR*"*"  (  Y,  Z  €  )  be  the 


eigenvectors  or  generalized  eigenvectors  of  Ah  corresponding  to  the  eigen¬ 
values  Xlt  Xg . X*.  Then— 

(i)  the  matrix  YrZ  is  symmetric, 

(ii)  if  y1  exists,  then  X  -  ZY~l  is  the  solution  of  (ARE)  such  that  the 

matrix  (£  -  WX)  has  the  eigenvalues  Xj,  Xj . X*. 


[Proof] : 

The  proof  was  first  given  by  Potter  (1966)  and  generalized  by  Martensson 
(1971). 


Lemma  2 : 
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There  exists  at  most  one  stabilizing  solution  to  (ARE). 
[Proof] : 

Let  X\  and  Xt  both  bet  the  stabilizing  solutions  to  (ARE), 


STXx  *  XxE  -  XxWXx  +  Q  *  0  (1) 

ErX 2  +  XaE  -  XtVXa  +  Q  *  0  (2) 

subtracting  (2)  from  (1)  yields 

et(xx  -  Xt)  +  (^!  -  Xg)r  -  +  XgWX2  *  o 

•which  may  be  rewritten  as 

(E  -  WXx)T(Xx  -  Xj)  +  (X,  -  Xa)(E  -  WX.j)  =  0  (3) 

Since  both  Xx  and  X2  are  stabilizing  solutions,  we  have 

Re[\(X  -  WXt)]  <0  Vie  i,2 . n 

and 

Tle[\j(E  -  WX2)]  <0  V  j'  e  1,2, . n 

From  Property  1,  we  conclude  that  (3)  has  a  unique  solution 

(Xj  -  Xf)  «  0  .  or  Xt  *  X*. 


Q.E.D. 

Now,  we  are  going  to  state  the  main  theorem  of  this  section  which  gives 
the  necessary  and  sufficient  conditions  for  the  existence  of  a  unique  stabiliz¬ 
ing  solution  of  (ARE).  Without  loss  of  generality,  we  will  assume  that 
We  GGt. 


Theorem  1 : 

The  stabilizability  of  ( E,G )  and  Re[X$(Aff)]  i*  0  (  v  i  =  1,2 . 2n)  is 

necessary  as  well  as  sufficient  for  the  existence  of  a  unique  stabilizing  solu- 
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tion  of  (ARB). 


[Proof] : 

( Sufficiency  ) : 

Suppose  (£\ff)  is  stabllizable  and  Re[MA#)]  *  0  for  all  t.  Let  the 


columns  of  the  matrix 


Yi 

rJ 


{Y\,Yi  e  Enx'n)  be  the  eigenvectors  or  general' 


teed  eigenvectors  corresponding  to  n  eigenvalues  -with  negative  real  parts 
and  J  be  the  corresponding  (real)  Jordan  block,  Le.. 


,r«. 

Yi 

J 


or 


EY ,  -  GGtYz*  YiJ  (4) 

-QY^EtY2^YzJ  (5) 


"We  will  first  prove  that  the  matrix  3^  is  nonsingular.  Suppose  Yx  is  singular 
and  %x&{Yi)  denotes  the  null  space  of  y».  Then  v  v  e  Ker^i),  we  have 
Y\V  =  0  and 

vtYzt  x  (4)  x  v  : 

~  vtyzteYiv  -  vryjGc!’y1!v  =  vTYtTriJlu 

Since  YXTYZ  is  symmetric  (from  Lemma  2-1). 

vTYiTY1Jv  svrYlrYiJv  =  {Yxv)TYzJv  s  0 
-»  GtYiu  -  0 

(4)  x  v  : 

■»  EYxv  -  GGtYzv  = 

->  rj/v  =  o 

i/v  £  KeriTi) 


(6) 
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It  is  clear  that  is  invariant  under  J  and  is  spanned  by  some  sub* 

set  of  the  (generalized)  eigenvectors  of  /.  Therefore,  3  eigen-pair  (  \,v)  ot  J 
such  that  Jv  *Xv  andv  C  Kar(Yi).  Then 

-QYxv  -ETY#i  =  YaJv 
-ETYzC  ~\Ytv 

Er{Y,v)  *  (-*)(}>)  (7) 

(-A)  is  an  eigenvalue  of  Et{ot  E) 

Furthermore,  we  know  that  Re(  -A)  >  0.  By  the  assumption  of  stabilize- 
bility,  we  conclude  that  ( -X)  must  be  controllable.  Thus,  from 
TheoremO. 2. 1.  l(iii)  (PHH  rank  tests) ,  the  matrix  [  -  E  \  G  ]  must  have 

full  rank  n. 

But,  from  (6)  and  (7),  we  have 

(Y&)T [  -£/  -  E  |  G]*0 

This  is  a  contradiction,  and  therefore,  Fj  must  be  nonsingular  and  .from 
Lemma  2-1,  we  know  X  -  Y%Yi~l  is  the  stabilizing  solution  Applying  Lemma 
2-2,  the  uniqueness  of  stabilizing  solution  is  guaranteed. 

( Necessity  )  : 

If  X  is  a  unique  stabilizing  solution,  then  ( E  -  GGTX)  is  asymptotically 
stable.  This  implies  that  (E,G)  mu3t  be  stabilizable.  because  of  the  sym¬ 
metry  of  MAh)  along  the  imaginary  euds.  we  conclude  that  #e[A<(Aff)]  *  0  for 
i  —  1,2, . ,2n. 


Q.E.D. 
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Corollary  1 : 

If  Q  ae  0,  then  the  stabilizing  solution  X  St  0. 

[Proof] : 

EtX  +  XE  -  XVX  +  ?*0 
(E  ~  VX)TX  +  X(E  -  VX)*  -(XVX  +  Q) 

(E  -  BT)  is  asympotically  stable  since  X  is  the  stabilizing  solution.  Thus, 
•the  solution  of  the  above  Lyapunov  equation  can  be  -written  as 

X  x  J  'l*‘K)T*(XVX  +  Q)  di 

o 

Since  the  matrix  (XVX  +  5)  •*  0.  X^  0  is  concluded. 

Q.E.D. 

Remark: 

The  proof  of  Theorem  1  -was  first  given  by  Kucera  in  1972  -with  the 
assumption  (?  se  0.  In  fact,  this  assumption  is  not  necessary. 

Theorem  2 : 

If  Q  -  HtH  i  0  in  (ARE)  and  X  is  its  solution,  then  Ker(j0  C  Ker (H). 
[Proof] : 

Since  A*  is  a  solution  of  (ARE),  we  have 

ETX  +  XE  -  XVX  +  HtH  =  0  (8) 

t  -a  ..  —  V t  VA  /  VI.  _  a\  'T*U  — 

LfZ.  V,  Ur  ^  AVUA  \SL  )  \yxwr  —  V/.  iUCU 

urx(B)x«; 

- utAtXu  *  utXAu  -  uTXWXu  *  urHrHu  =  0 


uTHTKu  =  0 
Hu  *  0 
u  C  K»r(H) 


Hence,  we  conclude  that  Ker(X)  c  Ker (H) . 
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2.3.4  Inner-Outer  and  Spectral  Factorization  : 


In  this  section,  the  special  form  of  coprims  factorizations  required  to 
reduce  the  general  H ,  optimal  control  problem  to  a  best  approximation 
problem  will  be  developed.  In  particular,  explicit  realizations  are  given  for 
coprime  factorizations  G  =  NM~l  with  inner  numerator  N  (Theorem  1)  and 
inner  denominator  M  (Theorem  3);  and  for  the  complementary  inner  factor 
TVjl  which  completes  the  inner  numerator  to  make  |a7  /fjj  square  and  inner 

(Theorem  2).  The  theorems  will  be  stated  for  right  coprime  factorizations 
(rc/)  with  the  duals  for  ic/'s  following  just  as  for  the  general  case  of 
coprime  factorization  developed  earlier. 


For  the  following  theorems,  it  is  assumed  that  G  - 


d! 

B 

Ci 

D 

e  1%"*  and 


the  realization  is  minimal.  We  will  denote  by  R^-  (i?5=0)  the  symmetric 
matrix  such  that  =  R  and  use  "D£  for  any  orthogonal  complement  of 
D  so  that  | DR~*  .Djj  (with  R-D'D)  is  square  and  orthogonal. 


Recall  from  Corollary  0.3.3. 1  that  N  - 


A_ 

\s\ 

C 

is  inner  if  and  only  if 


i)  S'X  +  d't  =  0  (1) 

ii)  t)'t)  =  /  (2) 


where  the  observability  gramian  X  solves 

A'X  +  XA  +  CC  =0 


(3) 


From  Lemma  2.4.1  a  stabilizing  state  feedback  F  yields  ref  G-NM~l 


where 


2 


U 

N 


A+BF 

BZ 

= 

F 

Z 

C+DF 

DZ 

(4) 


and  Z  can  be  any  nonsingular  matrix.  To  obtain  a  ref  with  N  inner,  we  aim* 
ply  need  to  use  equations  (l)-{4)  to  solve  for  F  and  Z.  This  yields  the  follow¬ 
ing  theorem: 


Theorem  1 : 

Assume  Jim,  Then,  there  exists  a  ref  G  -  NM~l  with  N  inner  if  and 
only  if  >  0  on  the  j  u-axis,  including  at  This  factorization  is  unique  up 
to  a  constant  unitary  multiple. 

A  particular  realization  for  the  factorization  is 


A+BF 

BR~* 

F 

R-* 

C+DF 

DR' * 

e 


RH lm 


where 


R  =  D'D  >  0 

F  -  -R~\B'X  +  D'C)  (6) 

and 


lA-BR-'D'C  -BR-'B1 

Xstac  [-CDiPiC  -{A-BR-'D'Cy  [  ) 

[Proof] : 

(only  if) : 

Suppose  G  -  Mr1  is  a  ref  and  NaN-I.  Then 
G*G  =  (. ~l)  =  >  0  on  the  ju~uds  since  M  e  RH~ 

(if): 


The  if  part  will  be  proven  by  showing  that  (l)-(4)  lead  directly  to  the 
above  realization  of  the  ref  of  G  with  inner  numerator.  That  G  -  NM~X  is  an 
ref  follows  immediately  from  (4)  once  it  is  established  that  F  is  a  stabilizing 
state  feedback.  Using  the  notation 

N  s  |  I 

we  will  use  (l)-(3)  to  get  N  inner  and  A+BF  stable.  From  (2)  we  have  that 
Z-R~^U  where  R~D'D> 0  and  U  is  any  orthogonal  matrix.  Take  U-I.  Equa¬ 
tion  (1)  implies  that 

R’*B’X  +  R-*D’(C+DF)  - 0 

so  solving  for  F  yields 


A+BF 

BZ 

C+DF 

DZ 

(8) 


F  S  -PrKB'X  *  D;C )  (9) 

Then  equation  (3)  yields 

0  «  A'X  *  XA1  -  t't 
=  {A+BF)'X  +  X(A+BF)  +  (C+DF)‘(C+DF) 

-  (A-BR-lD'C-BR~lB'X)'X  *  X(A-BR'lD'C~BR~lB‘X) 

*  (C-DR-lB'X-DR'‘lD,cy(C-DR-lB'X-DR~lD,C) 

=  (A-BR~lD‘CyX  +  X(A~BR~lD'C)  -  XBR~xB'X  +  CD^pfC  (10) 

since  DjP£  =  I -DR~XD‘.  Thus  Af=Ric  [An],  where 

\a~BR-1D’C  - BR-XB *  ,  . 

**  =  [-C'DjPiC  -( A-BR-XD'C )’] 


That  „3f=Hic  (4^)  exists  such  that  A+BF  ir  stable  folio1*?  from  Theorem 


2.1  as  follows.  Let 
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r  i 

CC  C'D 

(sI-A)~'B 

=  |B'(-s/-j4')m  I] 

D'C  D'D 

/ 

(12) 


That  T  =  G*G  >  0  is  true  by  assumption.  To  satisfy  Theorem  2.1  we  must 
have  (A  ,  B)  stabilizable  and  ( P  ,  A)  detectable,  but  this  is  immediate  since 
the  realization  for  G  was  assumed  minimal.  Thus,  Theorem  2.1  ensures  that 
X  -  Ric  (Ah)  exists  such  that  A+BF  is  stable. 


The  uniqueness  of  the  factorization  follows  from  coprimeness  and  N 
inner.  Suppose  that  G  -  =  NzMt1  are  two  right  coprime  factorize* 

tions  and  that  both  numerators  are  inner.  By  coprimeness,  these  two  factor¬ 
izations  .are  unique  up  to  a  right  multiple  which  is  a  unit  in  That  is. 


there  exists  a  unit  9  G  RH T*"1,  such  that 


9  = 


Clearly,  9  is  inner 


since  9*9  =  9*N*Ni9  =  N$Nz  *  I.  The  only  inner  units  in  RHm  are  constant 
matrices,  and  thus  the  desired  uniqueness  property  is  established.  Note  that 
the  nonuniqueness  is  contained  entirely  in  the  choice  of  a  particular  square 
root  of  R. 


Q.E.D. 


In  a  similar  manner  equations  (l)-(3)  can  be  used  to  obtain  the  comple¬ 
mentary  inner  factor  (CIF)  in  the  following  theorem. 

Theorem  2 : 


If  p  >  m  in  Theorem  1,  then  there  exists  a  CIF  Ni  G  RHJ,x^3~m^  such 
that  the  matrix  |jV  jVjj  is  square  and  inner.  A  particular  realization  is 


N  i  = 


U+BF  -FC'Dl 
C+DF  Dl  j  wi 


here  X  and  F  are  from  Theorem  1  and  X7  is  the 
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pseudo-inverse  of  X . 

[Proof]: 

The  proof  consists  of  verifying  directly  that  j//  //jj  is  inner  using  the 

above  realization  for  A^and  the  realization  for  N  from  Theorem  1.  Using  the 
notation 


N  tfil  s 


JJ 

B 

A+BF 

BR~*  -FCD\ 

ej 

7 

** 

C+DF 

DR' *  Dl 

(13) 


and  the  fact  that  Ker^)  C  Kbt(D£C)  (Theorem  3.2),  equations  (l)-(3)  follow 
immediately.  Thus  |/V  jVjj  is  inner. 


Theorem  3 : 


There  exists  a  ref  G  -  such  that  M  C  RH?™*  is  inner  if  and  only 
if  G  has  no  poles  on  the  j  oaxis.  A  particular  realization  is 


where 


A+BF 

B 

s 

F 

I 

,C+DF 

D 

(H) 


F=  -B-X 
and 


(15) 


X  =  Ric 


A  -BB' 
0  -A' 


i  0 


(16) 


[Proof] : 

The  proof  is  essentially  the  same  as  for  Theorem  1.  The  details  are 
straightforward  and  are  omitted. 
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Remark*  : 

(1)  The  minimality  condition  in  Theorem  3  can  be  weakened  to  ( A,B ) 
stabiiizable  and  A  has  no  eigenvalues  on  the  jotxis  and  the 
theorem  still  holds. 

(2)  If  G  €  RHmr>an  in  Theorem  1.  then  U  is  a  unit  in  RHm  and  Af“l  is 

"outer".  In  this  case,  G  ~  is  called  "inner-outer  factoriza¬ 

tion"  (IOF). 

(3)  Dual  results  for  all  factorizations  can  be  obtained  when  p  £  m.  In 

# 

these  factorizations,  output  injection  using  the  dual  Riccati  solution 
replaces  state  feedback  to  obtain  corresponding  left  factorizations. 


In  the  following  theorem,  we  may  assume  that  G(s)  is  stable  without  loss 
of  generality.  Any  G  e  RLm  may  be  factored  using  the  dual  of  Theorem  3  to 
obtain  a  stable  numerator  N  such  that  N*N  =  G*G. 

Theorem  4 :  (Spectral  Factorization) 

Assume  G(s)  €  and  y  >  ||C/(s)|U  Then,  there  exists  a 

M  e  RHmmMn  with  stable  inverse  such  that  MaIJ  -  y2/  -  GnG  with 


A 

B 

-R*KC 

R* 

where 

R  =7*/  “  D'D  >  0 
Kc  =  — R~\B‘X  -  D'C) 


X  S  KLc 


A+BR-lD'C 

C(I+DR-lD')C 


-BR-'B' 

(A+BR-'D’Cy 
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[Proof] : 


Let 


f  * 

- CC  -CD 

\sI-A)~xB 

-D'C  R 

1 

vhftre  R  *7*1 -D'D. 


Since  7>||C?!U  T(jt>)> 0.  The  minimality  of  the  realization  of  G(s) 

guarantees  that  (A  ,  S)  is  controllable  and  (-CC ,  j4)  is  observable.  Thus, 
from  Corollary  2.1,  there  exists  M(s)  €  Bp**7*  such  that  T  «  limM  and  a  par* 
ticular  realization  is 


A 

B 

Ji  * 

-R*Kt 

& 

where 


Kc  *  -R~l{B'X~D’C) 


and 


X  *  Ric 


A+BR-'D'C  -BR~lB' 
C(I+DR~lD')C  -(A+BR~lD'C)\ 


Since  G  is  stable,  we  conclude  that  U  e  RHmm1tm. 


Q.E.D. 
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2.3.5  Parametrizing  the  Optimal  Controller  and  the  Ha  solution : 


This  section  combines  the  results  of  Youla’s  parametrization  and  the 
coprime  factorization  to  parameterize  all  stabilizing  controllers  in  a  way  that 
is  convenient  for  solving  optimal  Lt  and  Lm  control  problems, 
let 


■fc 


Pi  i 

Pu 


A 

Bi  Bz 

- 

Ci 

Du  D\z 

Cz 

D\  i  Dzz 

Suppose  that  neither  P&  nor  P2\  bas  transmission  zeros  on  the  jcj-axis 
(including  *»  )  and.  without  loss  of  generality,  DuT *  /  and  DziDzi7  s  I. 
Under  these  assumptions,  let  (X>aa)j,  and  &L~  (Dzi)±  that  is.  |z?12 


h>r  'si 


and  li/21*  ,?/j_* 
and  H  given  as  follows  : 

F  =  ~U>i2 rCi  -  BzTX) 
\a -BzDizr  Ci 


are  orthogonal  matrices.  Them  factor  P  as  before  with  F 


x  =  mc 


-BzBzt 

-CSDgfCi  -(A-BzDjCiY 


and 


H  -  ~(BiDziT  -  YCZ7) 

\{A-BiDzSCtf  -CzTCz 
r=RlCl  -BiDlDgJ  -(A-BiJDzSCz) j 
Then.  Nu'Niz  =  /  and  P'zi^n9  -  I.  Also,  let  //j_and  Afj.be  CIF’s  30  that 

\a+b2f\Bz  -rcjD J 


\n„  n}  = 

l  "  *3 


c  ,  n  cl  n 


n.  I 
"i.  I 


A^HCz  \B\i-  HDz\ 

*  Cl 

= 

C2  1  D2 1 

**>.  j  <■*» 

2 


Letting 


T  * 

reduces 


(HV)n  N „ 

N21  0 


a+b2f 

-B2F 

Bx 

B* 

0 

A  +  HCt 

B  i  +  HD2x 

0 

Ci  +  DWF 

tDhF 

Dxx 

Dvz 

0 

Cz 

Bsx 

o. 

min 

XtX, 


I IL  | Ti(P,K)  stable 


to 


A 


a  min 

Q  C  RHm 


\\Ti(T.Q) 

(iVV)n  *  NuQNai 
the  optimal  K  may  be  recovered  from  Q 


Because  both  the  ||  •  ||g  and  ||  •  j|.  norms  are  unitary  invariant,  an  alter¬ 
native  expression  is  possible.  For  any  Q  e  RHa  (  a  s  2,  w  have 

j  (NV)u  +  RizQNzi 

ii 


=  i  [^12  ^j]“  ({ NV )u  *  NltQN2x) 


1 

J 

14 

\Vu*(N''?)nNvm  +  Q  Njtm(NV)uN£ 

NfiNViuNtS  N£(NV)UN£ 


Cjj 

On  . 


where  G  ~ 


G»  6*  *  [tffj  (W»  K’  ^ 


n.a  a  =  2  case  is  particularly  simple.  Since 


3 


Gll  +  Q  ^12 

' 

-  _!|2 

0  Gxz 

2  ' 

^21  ^22 

s 

p  , 

,1 

6n  +  £||a  + 

[Ggj  Gaaj 

2  , 

the  optimal  Q  is  seen  immediately  to  be 


«p* 


Unfortunately,  the  as*  case  is  more  complicated  and  will  be  investigated 
in  the  next  section.  To  obtain  an  explicit  expression  for  we  need  to 
compute  G. 


Note  that  (NV)n 


It  is  convenient  to  compute 


tfnmNu  ! 
.  *£  0. 


Vu 

r>* 


N- 


21 


seperately. 


□aim  1 : 


-{A*BzF)r 

(Cx-DxzF)TDlx  +  XBx 

-BtT 

D\ztDi\ 

vfcxr 

DfDn 

[Proof] : 


[*«" 

A7.» 

l*’X 


J  ‘ 


1^,!  = 


A  +  BZF 

0 

Bx 

-iCx  +  DxzF)T(Cx  +  DxzF) 

~(A+BzF)r 

-(Cx  +  DX2F)tDu 

n.-Ttr.  j.  n._r\ 

—  iC  \  1  •  — '  i«’  / 

B-t 

n.J’n.. 

Dfcx 

r  i 

-Dfcxr 

d£dxx 

conjugating  the  states  by 


I  0 
-X  / 


,  we  have 


At 


*-• 

W 

» 

r.f  i 

j4  +  ^  g/1 

0 

0 

Bx 

-(Cx  +  D„F)TDn-XB 

wi 

M = 

i?12r(C1  +  i?1s/,)+5Sr 

BuTBu 

0 

d£Du 

Since  Du 

•T(Ci  +  DitF)  +  B2rX  s  Di2TCi^BzTX^F  *  0  ,  the  claim  «  verified. 

Q.E.D. 


Claim  2 : 


•*{j4  +  HCz)7 

~  Ci?  y*B\D£ 

• 

s 

[Bx  +  HDzx  )T 
FY 

Dzit  Sf 

0  0 

[Proof] : 


Vil 

hi* 

conjugating  by 


i  y 
o  /'• 


A  +  tfC* 

-(Bx  +  HDzxfiBx  +  HDzx) 

(BrrHDzMhJ 

s,S/ 

0 

-{A+HCz)T 

Czr 

0 

-{Bx+HDzx)7 

ihxT 

/ 

0 

0 

0 

A  +  /fCj 

0 

0 

0 

0 

-U+^Cj)5, 

CbT 

-rs,£,r 

0 

i?!ir 

,  / 

FY 

0 

o 

which  verifies  the  claim. 


Putting  these  results  together  yields 


w. 


K 


Q.E.D. 


~(A  +  Bzr)r 

{C\  +  D\zF)T  D\\~  XB\ 

A+HCz)r 

-cj  rsiS tr 

a 

3zr 

DjOu 

* 

(£  i  +  HDsi)t 

Anr  d£ 

DfCJP 

DfDn 

- 

~FY 

0  0 

Note  that  this  is  the  cascade  of  two  systems  with  all  of  their  poles  in  C*  . 
Thus,  projection  onto  Hz  *  C  leaves  only  the  constant  term.  Therefore,  in 
the  case : 

Theorem : 


Q*  *Dx27Z-aDzit 


1 


2.4.2  Hank  el  Operators 

Let  G(s)  be  &  strictly  proper  transfer  matrix  which  is  analytic  in  Re  ssO, 
Le.,  totally  unstable.  The  Hanktl  operator  associated  with  G  will  be  denoted 
by  Tc  and  is  defined  as  follows.  Let 

P#  :  iaO'K)  -  (1) 

denote  the  orthogonal  projection.  Then 

?c  :  HdjR)  -*  HiijR)1 
Tcf  :=  PsiGf. 

There  is  a  corresponding  operator  in  the  time  domain.  Let  g  (t )  denote 
the  inverse  Laplace  transform  of  G(s)  and  let  Pg^  also  denote  the  orthogonal 

projection 

P„y-  Ia(R)  -  H2(TR)L  (2) 

(Context  distinguishes  the  two  projections  (1)  and  (2).)  The  time-domain 
Hankel  operator  is 

Tg  :  HS(1R)  -*  jy8(]R)i 

Tgf  ■-  PH±{g  +f  )• 

Thus 

u 

,  J^B(t-r)f{r)dr ,  f<0 
( Ttf)(t )  -  o ,  f  at  0. 


Since  the  Fourier  transform  establishes  the  isomorphisms 


ia(lR)  S  LzijTR) 
HtQR)  S  Hz(jTR) 
HgQR)L£  H&TR)L 


Z 


we  have  that 

INI  *  INI’ 


The  norm  of  T,  can  be  computed  by  state-space  methods  starting  from 
a  minimal  realization  of  G, 


G  * 


A 

B 

C 

o’ 

Let  Y  and  X  be  the  controllability  and  observability  gramians, 
o  o 


f  t*BB'*A,tdt  ,  X  :*  JtA'lCrCaMdt 


Lemma  1.  TjTt  and  YX  have  the  same  nonzero  eigenvalues.  In  particular 

IMI  *  M*- 

Proof.  For  this  proof  only,  drop  the  subscript  g  on  Tt  and  define 
.2.(0  :=  9  ("O'-  "Hie  adjoint  operator  of  T  is 

I*:  tfjOR)1-*  H,(R) 

r*’  '  ph±Isl*V- 

Let  be  an  eigenvelue  of  FT  anu  iet  /  e  Hg  be  a  corresponding  eigenvec¬ 
tor.  Then 

r T/  =  a3/.  0) 


Define 


3 


so  that 


Tf  *  ah  (+a) 

Vh.  s  gf.  (4b) 

(Vectors  (/  ,h)  satisfying  these  equations  form  a  Schmidt  pair  for  F.)  Since 

•  g(t)  s  C*MB,  t< 0 
&(t)  =  t>0 

we  get  from  (4)  that 

f  Ct^Bf  (r)dr  *  ah.(t),  f<0 
0 

0 

/ B't*MCh{T)d7  =  Gf(t),  t>0, 

% 

or,  equivalently, 

ftMv  s  ah(t),  t<0  (5a) 

B' t’rtw  *  c/(f),  f >0  (5b) 

% 

where 


v  :*  f t~ArBf  (r)dT 
o 

o 

w  :*  f  tA'TCh(r)dr. 

«*• 

Now  premultiply  (5a)  by  tA>t  C  and  integrate  from  —  •  to  0,  and  premultiply 
(5b)  by  A~~  B  and  integrate  from  0  to  •».  This  yields 


Xu  s  aw 


Yu i  =  av. 
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Finally,  we  get 

YXv  =  A,  (8) 

showing  that  ff*  is  an  eigenvalue  of  YX.  The  reverse  argument  leads  from  (8) 


to  (3). 


QED 


We  shall  be  concerned  with  approximating  G  by  a  stable  transfer  func¬ 
tion.  Le.,  one  analytic  in  Re  satO,  where  the  approximation  is  with  respect  to 
the  Lm  norm  Here  we  establish  only  that  the  distance  in  Lm(jTR)  from  G  to 
the  nearest  matrix  in  HJJTR)  equals  |f. 

Theorem  1. 


inf 


G-Q\  L :  Q  e  H~(jTR) 


*  llr' 


(8) 


and  tht  vnfvmxcm  is  achievsd. 

The  remainder  of  this  section  is  devoted  to  a  proof  of  Theorem  1:  only 
section  5  requires  some  of  the  material  to  be  presented  next 

The  inequality 


inf 


.  Ii~  j 

*  IM 


is  easy  to  establish.  Fix  Q  in  /?•(?' R)’  Then 


i  -  rll 

j  O’-vjj-  =  sup 


St  sup 


=  sup 


I  i/~  .  ■ !  .  —  ..  /  ■  .  1 1  .. 

; j2 :  /  ntjin) ,  (/ j 


\\PadG-Q)f\\  :/  e  HzijlR) ,  J  Ul 

II  -r 


! IHk1 


5 


*  INI- 

Take  the  inflmum  over  Q. 

It  is  convenient  at  this  stage  to  bring  in  Lz  and  Hz  with  respect  to  the 
unit  disk: 

T :  unit  circle. 

Lz(T,  Cmxn):  Hilbert  space  of  matrix-valued  functions  on  T,  with  inner  pro¬ 
duct 

<F.G> :=  ~^’trace|r(*^),C(e^)]d‘d. 

Hz(T,  Cmxn):  subspace  of  functions  F(z)  analytic  in  | z  |  <1  and  satisfying 

Lm(T,  C"1*"):  Banach  space  of  (essentially)  bounded  matrix-valued  func¬ 
tions,  with  norm 

\f\\.  :=  ess#sup 

Hm(T,  Cmxn):  subspace  of  functions  analytic  and  bounded  in  I  *  ;  <1. 


Map  the  right  half-plane  Re  s!eO  onto  the  unit  disk  |  z  |<1  via 


s 


s-1  _  1+z 

s  +  l  *  ‘  ~  l-z  • 


and  define 

G(z)  (10) 

♦ 

Since  G  in  analytic  in  Res  2:0,  including  the  point  at  <*,  G  is  analytic  in 
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\x  |asl,  i.e.,  G  e  Lm{T).  Associated  with  G  is  a  Hankel  operator.  Tj,  defined  as 
follows.  Again,  let  Pg^  denote  the  orthogonal  projection. 

PBi.  :  UT)  -  tfa(r>4 

Then 

Tz  :  H3(T)  -  H*T)L 

r zf  *  (id 


Lemma  2. 


Tc 


Proof.  Define  the  function 


i>(s)  *  V2/(s+l). 


The  relation  between  a  point  j  u  on  the  imaginary  axis  and  the  corresponding 
point  «**  on  the  unit  circle  is,  from  (9), 


This  yields 


di5  = 


\i>(ju)\*du. 


This  implies  that  the  mapping 

/  -  V/  :  Hz(T)  -  H&K). 


where  /(*)=/  (s)l 


is  an  isomorohism.  Similiarly. 


/  -  iff  ■  H3{T)L->  Hz(jTR)± 

is  an  isomorphism;  note  that  if  /  £  H3(  then  /  =0  at  z  =“,  so  that  /  =0  at 
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s=- 1.  and  hence  if/f  is  analytic  in  Re  s  <0. 

The  lemma  now  follows  from  the  commutative  diagram 

H2(T)  2  nt(T^ 

r« 

Hz(jTR)  -♦  JSTa(flR)l 

Q.E.D. 

There  is  a  matrix  representation  of  the  operator  Pg.  Let  the  power  series 
expansion  of  G  be 

G(x)  =  £  r‘Ct 

(Actually  the  sum  only  ranges  from  i*-«  to  i=0.)  In  (11)  let  the  power 
series  expansion  of  /  be 

/(*)  *  |z‘/t 

and  let  that  of  h.  :=  Ps^Gf  be 

A(z)  =  E 
■»«» 

Then  (11)  is  equivalent  to  the  equation 


h- 1 

G-i  G-t  ... 

u 

h.^2 

G-i  G-s  — 

fl 

h-  s 

G. a  C-4  ... 

f  z 

• 

•  • 

• 

a 

•  • 

• 

• 

•  • 

J 

• 

The  matix  in  (12)  is  the  familiar  Kankel  matrix  of  the  transfer  matrix  G{z ). 
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Proof  of  Theorem  1.  In  view  of  Lemma  2,  it  suffices  to  show  there  exists  a 
matrix  Q  in  Hm{T)  such  that 

|2-$||.  *  ||rj||.  (13) 

Let  the  power  series  expansion  of  Q  be 


<?(*)  «  £*<«. 

The  left  side  of  (13)  equals  the  norm  of  the  operator 
/  -*  (G-Q)f  :  HZ(T)  -*  Lt(T). 


The  matrix  representation  of  this  operator  is 


• 

t 

• 

• 

• 

Gj-9i 

£-<?. 

2-i 

2-i 

2-a 

2-s 

2-s 

• 

• 

• 

• 

• 

ft 

(H) 


The  idea  in  the  construction  of  a  Q  to  satisfy  (13)  is  to  select  Qt,Q |t...  in  turn 
to  minimize  the  norm  of  (14).  First,  choose  Qs  to  mirimize 


{  G-G  2. i  ...I  ! 

!  2_i  ...  i 

i 

•  • 

•  •  ... 


By  Parrott's  Theorem,  the  minimum  equals  the  norm  of  the  Kankel  matrix  in 
(12),  Le„  ,iTff!|.  Next,  choose  Q\  to  minimize 

if  •  i 
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Gx-Qi  G.-Q, 
G.-Q,  G-i 

£?-i  G- j 

•  ■ 

•  • 


•  • 


Again,  the  minimum  equals 


Continuing  in  this  way  gives  a  suitable  Q. 


Q.E.D. 


1 


2.4.4  Beat  Approximation 


The  transfer  matrix  R(s)  is  real-rational,  proper  and  enti-stable,  Le.. 
analytic  in  He  s< 0.  The  objective  is  to  find  Q(s),  real-rational,  proper,  end 
stable,  such  that  is  minimum,  Le.,  equals  ||r||.  The  constant  term 

of  R  can  be  absorbed  into  Q ,  so  we  can  assume  R  is  strictly  proper.  Further¬ 
more,  by  adding  rows  or  columns  of  zeros,  we  can  assume  R  is  square. 


Let 

R  * 
or  :* 

Q  * 


dj 

B 

PI 

Cl 

l 


(minimal) 


(to  be  found) 


Then 


G  :*  R  -  Q 


A  0 

B 

X 

o  2 

8 

c  -c 

-D 

A, 

B. 

c, 

A, 

So  Q  is  optimal  If  2  is  stable  and  |  |ffj  j_  a  a. 


Lemma  5.  Suppose  (3.P#)s.t. 


A,Pt  +  P,A,'  +  BtB%‘  =  0 

(5) 

C,P,  +  D,B,'  a  0 

(6) 

D,D,'  =  0*1 

(7) 

z 


P  it  0  -ukst*  P,  =  p  . 

Svppost  also  that  (A,P)  is  sitiblizabl*.  Then  Q  is  optimal. 
Proof.  Lemma  4  together  with  (5)  and  (6)  ■* 

G(s)G{-s)'  =  a2! 

*•>  !|c||.  a  ff.  Also,  (5)  ■* 

AP  +  M*  +  a  0. 

Thus  Lemma  3  ■*  A  stable. 


Recap:  objective  is  to  construct  A.B ,C \A,P ,£ ,3 ,P,  s.t. 


i) 

_  A  B 

E  9  pk 

(minimal) 

ii) 

A,Pt  +  P,A,'  -■  B,B,‘  a  0 

(5) 

iii) 

C,P,  *  D'Bt  *  0 

(6) 

iv) 

D,D,‘  a  (Pi 

(7) 

v) 

P*0 

(e) 

vi) 

(A.B)  stabilizable 

0) 

Construction. 

Step  1.  Find  a  balanced  realization  of  R: 


A 

B 

c 

0 

Thus  controllability  gramian 


(8) 


Q.E.D. 


=  observability  gramian 


3 


fflf  0 
0  2 


=:  P 


where  a  >  j  jz  1 1,  Le.,  r  »  multiplicity  of  a.  Partition  A.B.C  accordingly 

Un  ^1*1 


>4  = 


[Agj  An 


,  B  = 


B  i 
Bz 


.  c  =  c,  Ca 


Sttp  2.  Choose  3  s.t 


3B  i'  +  aCt  =  0 
33'  *  cr2/. 


(10) 

(ID 


Sttp  3.  Set 

B  *  -((Pl-X^-'CLBi  +  oCi'd)  (12) 

A  *  (-An* BtS')'  (13) 

C  *  CzZ  +  DBz  (14) 

P  a  2(ff2/-22)-1  (15) 


« 


P, 


- al  0  0 

0  -2  -I 

0  -/  P 


(16) 


Verification  of  (5)  -  (9). 

From  Step  1  we  have  that 

AP  +  PA'  *  BB'  -  0 
A'P  -  PA  -  CC  -  0. 

These  yield  the  following  equations: 


“  vA\\  +  B\B{  —  0 

(1?) 

—A\zZ  —  oAz  T  +  B\Bz  —  0 

(18) 

—An%  —  EAn  —  BzBz  —  0 

♦ 

(19) 

4 


mmoAii  “  oAjj  4  Cj'Cj  =  0 

(20) 

—.Ajj'll  ~  ffAjj  +  Cj'Cj  *  0 

(21) 

“.Ajg'L  “  EXjj  +  Cj'Cj  —  0. 

(22) 

To  see  that  Step  2  is  possible,  Le..  3  exists,  observe  from  (17)  and  (20) 

that 


B\B\  *  Cj'Cj. 

Hence  there  easts  a  unitary  matrix  U  such  that 


UB\  +  C\  *  0, 

Take  3  -  aU. 

To  verify  (5),  it  suffices  to  show  that  the  blocks  in  positions  (1,3),  (2,3), 
and  (3,3)  of  the  3x3  block  matrix  A%P»  +  P,A,  '  +  i?t  i?,  ’  are  rdl  zero. 


The  (l,3)-block  equals 


-Axa  +  BiS'  -  +  B1(B2'Z  +  c3Cd('Z*-o*I)"1 

*  +  (BiBi'Z  -  ffCj'Cj)(E*  -  a2/)”1 

=  “^12  + 

UijZ  ▼  ff/l8l')2  -  o(j4jj'E  +  <rAja)j(2a -a3/) 


from  (12) 
from  (10) 
from  (18);(21) 

from  (18):(2l) 


=  0. 


The  (2,3)-block  equals  zero  immediately  from  (13).  Finally,  the  (3,3)*block 
equals: 

AP  +  PA'  +  §§•  =  (-Aa  +  BBz’)P  +  P(-Azi  + B2B')  +  §§•  from  (13) 
=  -Aa'Zia2!-!*)-1  -  (o*/-S*)-*L4B  from  (12);(15) 
-(o2/-22)-l(25a  +  ffC,a'j5)5T(<72-Sar1 
-(o2/-2Si)-1252(5?  T  +  oS’CzXo2!-#)-' 
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+(<ja/-28)-l(2£a  +  aCz’£)(B31  + 

«  (ff2/-2®)”1(~esAila'2  +  -£S252'E 

o^EAja  +  SrfljsS*  +  ff^Cy  C^Xo8/— E8)”1 
•  0.  from  (19);(22). 


Next  is  the  verification  of  (8).  We  have 

C,Pt  +&,£,•  *  “35s’~C82^&,  -W'-C*-5^] 


sc 


p.  0. 


from  (10)  and  (14). 


And 


J5 S'  *  C2  +  «  0 

by  .substitution  from  (12).  (14).  and  (15). 

Equation  (7)  is  immediate  from  (11),  and  (8)  follows  from  the  definition 
(15)  of  P. 


It  remains  to  prove  (9).  Suppose  there  exists  A,  He  A5*Q,  such  that 

[a-w.*]' 


doesn’t  have  independent  rows,  i.e., 


‘[a -A/, 


B 


a  0 


for  some  s^O.  From  (13)  we  get 

*•(-■<*«’  +  &Bt')  =  As*,  s'jg  =  0, 


so  that 

f — r*  "7 ■» 

\  ***41—  * - 

This  implies  that  (-Azz)  is  unstable,  which  is  not  possible:  stability  of  («vi) 
implies  that  of  (-Aw). 


1 


2.4.5  General  Distance  Formula 


Consider  the  problem  of  minimizing 


R 13 

Rzi  R» 


where 


R  := 


Ru  R ia 
R 2i  Ra. 


is  real-rational,  strictly  proper,  and  analytic  in  Re  ssCO,  and  Q  is  required  to 
be  real-rational,  proper,  and  analytic  in  Re  siO.  This  section  contains  a  for¬ 
mula  for  the  minimum  h  w.ns  of  the  norm  of  a  certain  operator  F.  Note 
that  the  minimum  is  the  distance 


I RHm 

Mo 


dist  (R 

from  R  to  the  set  of  all  matrices  of  the  form 


0 

0  ) 


Q  0 
0  0 


L  Q  e  RHm. 


The  matrix  R  induces  a  multiplication  operator  on  Lz  @  Lz,  Le„ 


/ 

V 


Ru  R\t  f  i 
R a:  Rzz \y  z\ 


r  :  Hz  ®  Lz  **  Hz  ®  Lz 


\fi\  P*° 


o  / 


Ru  Rizlfi 
Rz i  Rzz[  f  z ' 


Define  F  via 


2 


Theorem  1 


dist(J? 


RH. 
'  0 


0 

0 


The  proof  of  this  theorem  is  a  straightforward  application  of  the 
theorem  of  Parrott  /  Davis  -  Kahen  -  Weinberger,  as  in  Section  2.  It  is  possi¬ 
ble  to  write  down  a  matrix  representation  of  I\  but  an  efficient  numerical 
procedure  for  computing  its  norm  has  not  yet  been  developed. 


STF  (M’S 


M  ores 


ONR/HONEYWELL  WORKSHOP 


Problem  Description 
and  Motivation 


Real  Control  Design  Problems 


^  a r«  ^  /*\  m*\ <»ty  a  ^ 

V'yMii^vii'ym.g  v.v  mur\c  u i<5  ^iuiiw 


"behave  well* 


Control  Designers  Need  Many  Tools 


Simulation 


Identification 
Methods 


Formal  X) 

Synthesis 

Problems 


Software 

Development 

Aids 

Laws  of  Physics 


r^^Classical  Feedback 
Theory 


Numerical 

Optimization 

Methods 


Formal  Synthesis  Problems 


•  Some  Element  of  Design  Abstracted 
as  a  Formal  Mathematics  Problem 

»  Solutions  Known  Mathematically  and 
Computable  Practically 

■  Properties  of  Solutions  Well-Understood 
and  Desirable 


-The  closer  the  abstracted  problem 
matches  the  real  design  task,  the 
more  powerful  it  is  as  a  design  tool 

—No  abstracted  problem  matches  the 
real  thing  perfectly 


Cur  Formal  Synthesis  Problem 


plant,  sensors  and  actuators 
abstacted  as  a  linear  time-invariant 
finite-dimensional  operator 


measurements 


computer  and  software 
abstracted  as  another  linear  time- 
invariant  finite-dimensional 
operator 


Design  Task 


Specify  K(y)to 


■  achieve  stability 

•  optimize  performance 

•  provide  robustness 


Stability 


many 

other 

outputs 


No  bounded  input  should  produce 
unbounded  outputs 

- — >  we  must  achieve  flctorasll  sfLaMfliiy 

(all  poles  in  the  left  hand  plane) 
(stabilizable/detectable  design  models 


Performance 


selected 

inputs 


selected 

responses 


Inputs  in  a 

Specified  should  produce 
Class 


Responses  with 
Specified  Optimal 
Properties 


Performance  Objectives  for 
A  Familiar  Example 


•  Input  Class 


Sample  Functions  of  a  Stationery  Random 
Process  with  Specified  Autocorrelation 


white 

noise 


linear  operator  - 
staple  and  minimum  phase 


v* 


k 


•  Response  Properties 

« 

Minimum  Mean-Square-Error 


selected 
responses  — 


min  E 


linear  operator  -  ste&le  end 
_  mlmmum  phase 


{  e(t)  eCt)  } 


Resulting  Design  Model 


Overall  Plant: 


**  «a 

e 

P  P 

m  c 

n 

-  P  (n,u)  = 

11  12 

Lj 

P  P 

L  21  22  J 

u 

n  « 

Closed  Loop  Responses: 


e 


+  P..II-KR 


12 


22 


J"1  KR 


21 


n 


r(P.K)  n 


Resulting  "  H2'Optimization  Problem 


I  flin.  E  {  e(t)Te(t)  } 
|  stabilizing  K 


ECe(t)  Te(t) }  =  E  {  [  F(P,K)  n](t)  [  F(P,K)  n](t)  } 


=  Tr  {  J  F(t)  F(t) T  dt  ] 
0 


'impulse  response  of  operator  F 
=  Tr[^r\  F(jw)F(jw)*  dw  } 


-oo 


F(s) 


frequency  response  of  operator  F 
2 


F(s) 


=  *]/ T r  J  F(]w)F(jw)  dw 


oo 


2-NORfl  on  the  Hardy  space 
of  stable  transfer  functions 


An  Equivalent  Deterministic 
Performance  Objective 


e(t) 


Input  Class 

m  specific  time  functions  in  L2  representable 
as  impulse  responses  of  linear  systems 


impulse  j 

into  j-th  _ 

channel 


linear  operator  - 
stable  and  min  phase 


•  Response  Properties 

Minimum  L2“Norm  of  m  corresponding 
weighted  responses 


selected 

responses 


linear  operator  -  stable  and 

minimum  phase 


min 


2 


eJ(t) 


2 

2 


e(t)Te(t)  dt 


2-NORM  on  the  space 
of  time  functions 


Alternate  Performance  Objectives 


•  Input  Class 


All  possible 
Injunctions  with 
bounded  norm 


II v(t)llj  si 


•  Response  Properties 

Minimum  L2Norm  of  the  largest  weighted 

response 

selected 
responses 


SUp 

V 


2 


min 


Resultirg  "  H^-Optimization  Problem 


— 

Min  S“P 

e(t) 

stabilizing  K  v 

2  | 

sup 

V 


sup 

V 


FCP.K)v(t) 


Cf  D 
» 


V  \ 
fi'S 


■L2 


operator  norm  induced  by  the 
2-norm  on  input/output  time  functions 


oo-Norm  on  the  space  of 
stable  transfer  functions 


More  General  Performance  Objectives  ? 


All  possible 

IL  -inputs  with 
P 

ilv(t)  iip<  1 


> 


V 

ttttnwtsiu 
it  !| 

iw  H _ ■ 

j  i 

j  G 

r 

ij  ji 

L J!\a/  i- 

e 

u 

7;  W  •  jj— * 

|i  '  \i 

tSiStWtiiUSt* 

i  !woir 

y 

i 

_ J 

\s 

l: _ 

K 

T - 

Minimum  IL  -noi 
of  the  largest 
response 


■  Trivia]  solutions  (optF=0)  whenever  r  >  p 

■  Little  engineering  precedent  for 

3  <  r  <  p  <  oo 

•  Bounded  signal  problems  (r=p=<») 

indirectly  covered  by  square-integrable 
problems 


For  5150  Operators  (Doyle/Gahberg) : 


F(t)||  <  2M  F 
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Me  Millan  Degree 


Robustness 


Stability  must  be  maintained  in 
presence  of  L 


Some  minimum  performance  level 
must  be  maintained  in  the  presence  of  L 


Some  Philosophy  About  Perturbations 


Nature's  Perturbations 

•  are  unknown, 

•  are  potentially  catastrophic  for 
any  control,  and 

'•  defy  mathematical  description  . 

Nevertheless,  we  must 

■  represent  them  by  mathematical 
models , 

« 

•  specify  maximum  levels  of 
severity  within  those  models ,  and 

•  desiqn  controls  to  work  successfully 
for  this  specification 


The  “Leap  of  Faith"  that  the  selected  representation 
will  protect  against  the  real  thing  remains  the  burden 
of  engineers,  not  of  mathematicians 


A  Model  of 

"Unstructured M  Perturbations 


L  (z’)  =  W  A  W  (z'J 

^-stable  minimum  phase 
f  LFDTI  weighting  operators 

any  LFDTI  operator  from  a  norm-bounded  set 

p  ={  A  A  stable,  IIa|I  <  *} 
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Transfer  Functions  in  a  Disk:  cr[A(jw)]  <5 


log  gain 


Design  Model 


Overall  Plant  P 


Robust  Stability 


System  S  : 


Theorem 

Given  :  i.  F(P.K)  stable 

ii.  A  <•  (P  [stable  with  ||A||^2  £  5] 

Then  systemS  is  stable  if  and  only  if  ||F||.  <  ~ 

o 


Resulting  Analysis  Test  for  Robust 
Stability: 


NOTE:  This  is  8 
non-conservative 
test  for  stability 
with  respect  to  the 
perturbation 
set  (p 


Other  Stability  Robustness  Tests 


Set  of  Perturbations 


Individual  scalar  perturbations 
acting  one  at  a  time 


<p  =  (a  I  A=diag  (  0,..0,  Aj,0,..0  ) 
1  *■  I  Aj  stable,  liAjllSS 


A  single  scalar  perturbation 
acting  simultaneously  in  each  loop 


p  =  [A  A=diag(  A, A,..  A  ) 
1d  ^  A  stable,  II A  ii  <  S 


•  A  single  multivariable  perturbations 
=  {A  |  A  stable,  II  A  it  <  5 


} 

} 

} 


Corresponding 
Robustness  Tes 

|  F..(jw)|<—  w  £ 
for  eacl 


p[F(jw)]<~  w  < 
o'  [  F(jw )]  <  y  w  < 


•  m  multivariable  perturbations 
acting  simultaneously 


x=  {a 


a 

thog 

Aj  €  ( P 


<,A,.A,,...A  J 


P  [  F(ju)]<  J  w  < 


Conservatism  ? 


»  Tests  are  conservative  (sufficient  but  not 
necessary)  whenever 


set  of  true 
perturbations 


C 


f  set  of  modeled 
perturbations 


•  Simple  norm-bounded  covering  sets 
are  often  conservative 

Examples 

-  Real  parameter  variations 

-  Deliberately  neglected  dynamics 

(time  delays,  high  freq  lags  ) 

-  etc 

■  Of  norm-bounded  sets,  induced  2-norm  ( IHloo  ) 
sets  are  least  conservative 

II MU.  <  H Mil. 
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For  5150  Operators : 


Robust  Performance 


p— If—— a 
—  £[<• 


Robust  Performance  Specification: 


for  all  A  C  <P 


Robust  Performance  Specs 
Viewed  as  Stability  Conditions 


K'&MmmmBmmmmammmmmmmmmmmammmmmmiixmmcxKwemmmiwmvMwimmx 


Theorem 


||FjFi(P.IC>.A]||Aa  <1 

iff  (i)  fJFjCp.kj.a]  is  stable,  and 

(2)  the  following  system  is  stable 
for  all  Aq  €iP 


o 

<] 

f- - 

— “> 

fJf^(P,K),A] 

(stability  robustness  theorem  used  backwards) 


Resulting  Condition  for 
Robust  Stability  ©cod  Robust  Performance 


Fu[F^(P,K),a]  stable  for  all  ZS  €  (P , 


and 


stable  for  all  Aq  6  (P 


....  jji|  ill 


stable  for  all  Aq,  S  tf3 


Stability  and  performance  robustness 
are  achieved  simultaneously  if  and  only 
if  our  feedback  system  is  stable  for  all 
perturbations  with  a  particular  2x2 
block-diagonal  structure 


Once  More  for  Emphasis  .  .  .  . 


•  Conditions  for  Robust  Stability  and  Robust 
Performance  are  Equivalent  in  the  induced 
2-Norm  (H  ) 

oo 


A(s) 

Stable  A  € 

F(s) 

Fts) 


OO 


J_ 

5 


»  Simultaneous  Stability  and  Performance 
are  guaranteed  by  a  Stability  Condition 


0 

0  A 
*•  •# 

Stable  A„ 

1  u 

A  €  P 

Fts) 

A  Conservative  Analysis  Test 
for  Robust  Stability  and  Performance 


System  S 


1 - 1 

f>  o 

ll _ 1 

/ _ 

r 

WJ 

_v _ 

Fj(P,K) 

2 

e 

J 

remains  stable  for  all  A0  .A  6  <P  ir 
a  l  Fj  (jw )  ]  <  -J  for  all  w  s  ~ 

(standard  stability  robustnesslheorem 
with  structure  ignored  ) 


E®DQ@©(FV©ftil^© 

cr [  Fg  ( jw ) ]  =  |f  (jw  )| 


Example:  = 


0  f 
0  0 


L  I.  ‘1  • 


lest  gives  staomty  only  for  sup] f  (jw  ) j 

w 

Actually,  stability  is  never  lost 


M)  -det 


1  -fA 
0  1 


*  0 


A  Tight  Stability  Test 
for  Structured  Perturbations 


■  Set  of  Structured  Perturbations 


x  =  {a 


a-SI5kia1.a2 . AJ  I 

Aj  €  <P  J 


A(s)  =  ^flA1(S)J...,Am(S)] 

O'  [ A  (jw)  ]  <  &  •  for  all  w  and  all 

J 


} 


•  A  New  Function  of  a  Matrix  : 

The  "Structured  Singular  Value  (5SY)~ 

1 

ptF(iw)]=  mjnf5  det[l-F(jw)A(jw)]=  0  j 

*  for  some  w  and  ^(6)  ^ 


Theorem 


Svstem  S 

1  - 

is  stable  for  all 


A  €  X(<5 )  iff  ji[F(jw)]  for  all  w<  oo 


Mu  Makes  a  Difference  ! 


Example 


A,  0 

Lo  *2 


cr [  F(jw)]  =|f  (jw)| 


Singular  Yalue  Test  gives  stability  for  s£f|f(jw)|< 

Structured  Singular  Yalue  Test : 

1 

=  rjjjf,  I"  5  I  det[I-F(jw)A(jw)]=  0  j 
*  I  for  some  w  and  A  6  X(6 )  * 

1 


min  (  5  det 

' i  -f(jw)  t~2 

Ia  | « r 

o 

1 

oo 


Stability  is  maintained  for  all 


Mu  is  a  Potent  Tooi  ! 


An  "  H ^-Optimization  Problem  " 


Min  sup  y  [F(P,K)(jw)]  | 

stobilizing  K  w  4 


•»«•«* . . 


'.V«VrtV«* 


•-Optimized  Robust  Stability  and 
Performance 

s“p  JJ  =  Popt  ==>  Stability  for  all 

A  €  X  (—  ) 

^opt 


Guaranteed  Performance  Level 


Other  "IHljj- Problems" 


Maximized  Performance  5ubject  to 
Stability  Robustness  Constraints 


„  UJ  s.u,p  [  F  ( P,K,6_,  6  =  1  ]  <  1  } 


j  stabilizing  K  v  0  '  W  "  *  - 0J 


■  Maximized  Stability  Robustness  Subject 
to  Performance  Constraints 


Max 

stabilizing 


(« 


sup 

w 


U[F(P,K,60=  1,6]  <  t  } 


Summary  of  Formal  Synthesis  Problems 


Class  of  Performance  Set  of  Synthesis 

Inputs  Objective  Perturbations  Problem 


White  Noise  E(eTe)  A  =  0 

Responses  A -0  j 

>  Min 

F(s) 

2 

l  -bounded  ,  . 

Signals  sup  !1.e(t)  "2  A=0  ] 

Min 

F(s) 

oo 

-  -  A  €  <P  J 

(  sup  CT[  F(jw)] } 

OJ 

H  -bounded  ^ 

Signals  sup  lle(t)  1!2  A** 

Min  | 

V 

F(s) 

1 

'll 

{  sup  \x{  F(jw)j  j 

CO 


At  X  J 


ONR/HONEYWELL  WORKSHOP 


Historical  Perspective 


A  Generic  Classical  Control  Problem 


All  SISO 


Perturbations 


controls 


Design  Objectives  Expressed  in  Terms 
of  Desired 


Loop  Shapes 


Small  outputs  in  response  to  disturbances 


ll+GWjw)]'1  Wd(jw) 
GK(j  w ) 


<  1 

>  |Wd(jw)|  GK>>  1 


•  Stability  in  the  face  of  perturbations 
j  GK[l+GK(iw)  V]  W,  (jw)  j  <  1 

|GK(iw)l  <  G.K<<1 


Desired  Loop  Shapes 


Multivariable  Generalization 
cf  Loop  Shaping  idea's 


\ 


Interpretation  of  Loop  Shaping 
In  Terms  of  General  Design  Model 


SSV  for  the  Classical  Loop 


■  Definition 

1  * 

=  minrs|  detl I- Fljw)A(jw)]  =  0  j 
t  I  for  some  w  and  A  €  X(5 )  ' 

A  r  a,  o  i 

with  A=[  0  Ad. 

•  Some  Calculations 

1-f.a  =  i-  [GKti+eKr;w'i  ca>  a<o 

*  |  (l+GK)"  Wd 

is  singular  iff 

1  -  A,GK(I+GK)''  W,  -  Ad(!+GK)"'  W„ 
is  singular 


<  8 (gko+gk)"1  Wj]  +  er [ti+GK)-1  wd ] 


Loop  Shapes  which  satisfy 
classical  objectives  tend  to 
minimize  ji 


Limitations  of  Loop  Shaping 


|  All  design  objectives  must  jj 
j  be  reflected  to  one  point  in  |j 
||  the  loop  j 


Loop  Shaping 
works  for  this 
■multivariable 
loop  . . . 


What  Goes  Wrong  ? 


i—  l+L  - 

G 

-o — 

-K 

..... 

w 


I*'*' 

1 

*1- 

*d 

| 

Fj(P,IO  = 

KGd+KGr'w  Kd+GKr'w. 

1  0 

Gd+KGl'V,  (l+GK)"'w 

1  Q 


U 


-i-KH 


-K 


m 

W[ 

I 

•  • 

[sd+KGj'Vj  d+GKr'w,  ] 


I-F^A  is  nonsingular  iff 


I  - G  (l+KG)”  V  A  K  -  (l+GK)“ '\V\A„  is  nonsingular 

A  II  a  d 

K^K 


r1< 


y E  < 

I  Classical  objectives  do  not  minimize 
jj  unless  condition  of  K  remains  small 


CTK  — 

^jcrtKGd+KG)'1  W,]  +  cr[(l+GK) 


wdi; 


An  Aside:  The  Scaling  implicit  in  ]1 


Fg(P,K)  = 


KGd+KGrV  K(l+GK)’1Wd 
Gd+KG)"1^  (l+GK)"1Wd 


K 

I 


(l+GK)"1  [GW,  Wd  ] 


•  To  compute  p,  find  the  minimizing  D 
for  5  [ DFZ  D“ 1  ] 


D^D”1 


(1+lGdK)’1 

d 


D  = 


"dl  0" 
0  I 


%  Equivalent  to  a  change  of  "units’* 
in  the  problem 


Decoupled  Flight  Path  /  Attitude  Control 


Decoupled  Flight  Path  /  Attitude  Control 
In  Our  General  Framework 


\ 


Selected  Weighting  Functions 


An  Aside  :  Weights  are  the  "Knobs" 

of  the  Formal  Synthesis  Problem 


Selected  Desired  Loop  Shape 


5j(KG) 


implications 


mmrmxmiM9wntm*ruLmumit 


•  The  Loop-Shaping  Design  is  Robustly  Stabie 


W.K6  (l+KG) 


OO 


1 


•  Performance  is  Satisfactory  for  the 
Nominal  Plant 

||  W)(l+GK>'  ||  <  ] 

•  But  there  exists  a  A,  with  | Ja, |  <  1/15 

such  that 

||w2Fu[FE(P,K).A1]|L  >  15 

®  The  design  does  not  have  robust 

•  m  * _ 

per  1  US  HI£25lC8  I 


ONFi, 'HONEYWELL  WORKSHOP 


Overview  of  Synthesis 
Solutions 


0[Wi(f+KG)-'KG) 


Summary  of  Formal  Synthesis  Problems 


Class  of  Performance  Set  of  Synthesis 

Inputs  Objective  Perturbations  Problem 


White  Noise 

Impulse 

Responses 

E(eTe) 
Illei(t)  ll2 

r  \ 

O  O 

II  II 

<  < 

0 

Min  ||  F(s)  | 

> 

{  J  T  r[  F  *  F  I  i 

2 

J/2 

dwj 

Unbounded 

Signals 

sup  lle(t)  il2 

A  =  o  1 

i 

Min  | 

F(s)  | 

oo 

— 

— 

A  €  <P 

✓ 

■ 

{  sup  CT[  F(jw)j} 
to 

Unbounded 

Signals 

sup  lle(t)  ll2 

Min  | 

s 

r 

F(s) 

p 

• — 

— 

A«  X 

(  sup  F(jw 
to 

■n } 

Status  of  Synthesis  Solutions 


■  The  ^"Problem 

Solutions  completely  known 
(Parallel  Wiener-Hopf  /  LQG  Solution  ) 

■  The  IHl^Problem 

Solutions  available  arbitrarily  close 
to  optimal 

»  The  Ini  -Problem 

Approximate  Solutions  available  through 
Iteration 


The  W  -Solution  Process  a=2,oo 

oc 


Min 

F.(P,K) 

I  Stabilizing  K 

z 

OC  : 

a  p  z 


K 


Parametrlzeuion 


Min 

Qtm 


if 
Hi* 
l|f; 

if 

:|il 

jj>i 

fi-iiiijijij1'- 

";i||i!i  ‘ 


F  ■ 


r”1 


P. .+  R  ,K(I-P  K)  P, 
11  12  22  21 


WT2, 


V2V  V*.“ 


jiji 


Unitary  Invariance  $ 


Min 

Q€RW 


oc 


*VQ  R12 
R21  R22 


Projection 


oc 

Dilation 


Min 

Q€RM. 


Ri  t+^  'I  2 


Min 

Q€RM 


lie +8  || 


oo 


00 


°  =  -pH2(r„) 


:s;i  ft.. 


G  =  Best  Approx 

..of  L  G 


fjfii 


Re-Parametrlzation .  P" 


K 


opt 


Parametrization  Highlights 


e 


t 


T(s) 


e 

Residua] 


T22(s)  =  (residual)  /  u  =  0  (  separation  theorem  ) 


F^a,Q) 


Qt  l-T  0) 
22 


Tn+T12QT2, 


Parametrization  Highlights 


M 


•  Youla  Parametrization  of 
all  Stabilizing  Controllers 

K(s)  =  (  U  +  MQ  )  (  V  +  NQ  )*’ 

0  0 

where  P22=  NM  KQ=  UQV0 1  stabilizes  ,  Q  stable* 

z 

•  Observer-Based  Compensator  (OBC) 
with  added  stable  dynamics 


Exploiting  Unitary  Invariance 


21 

JL 

21 


VVT2,  * 

= 

[T.*Tt]*[T, 

R,,+Q  r,2 

i 

R  R  ^ 

»  *  21  '  22 

/ 

THE  PRICE  : 

unstable  functions 

* 


Projections 


Rn+Q  R,2 

R2I  R22 


oo 


Tr(R  +Q).  (R  +Q)+  Tr(R*  R  )  +  Tr(R*  R  )  +  Tr(R*  R 
11  12  12  21  21  22  22 

-oo 


It  is  sufficient  to  minimize  Hr  +q 


i  i 


turns  out  to 
be  entirety  unstable 


must  be 
stable 


Q  =  0 
opt 


)  du> 


•  ■ 


<REE : 


Note  5 


interconnection 

Matrix 
Coniro  Her 

topics 

i .  we  U  -poseJness 
Z.  internal  stability 

3.  parametri nation  of  K 

4.  state  -  space  realiesfion 

5\  closed -loop  transfer  matrix 
key  references 

I.  Yam  I  a  .  -Tafer-  .  SonOt  l  etsrn/y  '*7  Ct 

"  — j  -  —  #  -  ^  -  •  ■ 

z,  Oesoer  y  LiUy  Murray }  Saelcs  1  &o 

♦ 

3.  Me/f,  Jeco  bs&n  y  Salas  l*4 


oTAglLl  ZATION 


p  drdtnC'f  ri2d{  ion  o{  K 


(.  wt  ll- po:cJ»ess 

2.  internal  staliliiy 

3.  o9tam.e{rii-atiom.  of  K 

4.  $  late -space  realieafio* 

s.  closed  ‘loop  transfer  m affix 


Assume  P ,  K  proper 


wow*  A dirioMS  exogenous  inputs 
Jef'n  system  is  well-posed  tf 

exists  X  is  proper 


4rahst<r  matrix 


suffices  to  consider 


u 


V  (s  proper 


Cone  (us  ion 

well- pose  J  ness 

4==*  I  -  K  W  Pu  W  invertible 


4=r=rr» 


4r-  fr 


I  -  P  {«©)  K  («*)  » t  I* 


r 


-K(*0 


invert  it  If 


well- pose  J  ness 


assume*/  heredfi  tr 


1.  well-/ toseJmeet 
Z.  internal  STali liiy 

3.  0jr»M^frii9lie*t  of  K 

4.  state  <■  space  realisation 

£  claseJ-!oop  transfer  matrix 
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C 

i 
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C  ( Si  -  aV'b 
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a 


def’*  infernal  slsLihiy  mignS 

*«),*«)—*  o  V  i.c. 

cLdrzticriddiion 


minimal  r£a/»?a/io*S 
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rA 

K  I 

p= 

c, 

P. 

Pa 

Cx 

m 

P« 

Pit. 

K  _ 

A 

A 

A  1 
8 
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Y*Ca**DMtt 
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[» v  *  *'ai 

KM  *  D 


elit*i*»4e  «,  y  ,  ,  „  < 

/  exins  by  wtfu- 


r  A ' 
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m 

i  -n 

0  c 

* 

a 

LYJ 

-K-  ij  ' 

« 

r* 

*• 

O 

i _ 

* 

*  « 

i**(6r*dl  sisLiliiy 


ues 


existence 


m 

A 

8, 

*il 

p= 

c, 

P„ 

Oa 

Ia 

K 

on. 

\tmn%*  3  proper  K  achieving 

inter**!  stability  *  * 

(At6%)  *idlili?*Ue 

(unstable  modes  ot  P 

CO*tr*ll*L(i  Prom  t* 

i  (Ct^A)  cKetectoUe 

(«*sUle  modes  *{  P 

obserudLle  *i  y) 

...  s$s***eJl  here  dfter 


sleklwLle  2 
deled  « l  (e 


leu*** a  left  lytle**  ider**tiy  eivKe 
«  rl9M 


**•  suffices  io  slalili&e  £ 

«i 


1.  well-feseJnesf 

2.  i *\ier*9  l  sfsLilify 

3.  psrmeiritdion  of  K 

4  .slaie- sp»ee  reakiaho* 
5\  closed  -  loop  IraAf/Vr 


*f  2cfcriE3^i0H$  over  R 


*f  6  R  H 


real  **  rational  3 
proper  p 

analytic  in  tie  $  &0 


.•.  extend  to  vectors  &  weir  ice  % 


de$'r\  n3  N  e  ar#  r » *|  W  -  C<rp  r  i  m  € 

[n]  '*  (eft  •  invertible  in  RH^ 


equivalently 
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